Methods for determining relative phase differences in entangled quantum states

ABSTRACT

Various embodiments of the present invention are directed to methods for determining a phase shift acquired by an entangled N-qubit system represented by a NOON state. In one embodiment, a probe electromagnetic field is coupled with each qubit system. The phase shift acquired by the qubit systems is transferred to the probe electromagnetic field by transforming each qubit-system state into a linear superposition of qubit basis states. An intensity measurement is performed on the probe electromagnetic field in order to obtain a corresponding measurement result. A counter associated with a measurement-result interval is incremented, based on the measurement result falling within the measurement-result interval. A frequency distribution is produced by normalizing the counter associated with each measurement-result interval for a number of trials. The phase shift is determined by fitting a probability distribution associated with the probe electromagnetic field to the frequency distribution as a function of the phase shift.

TECHNICAL FIELD

The present invention relates to quantum computing, quantum information, and quantum metrology, and, in particular, to methods for determining a phase shift acquired by entangled, quantum systems.

BACKGROUND OF THE INVENTION

Advancements in quantum physics have facilitated understanding and enhancement of existing technologies. For example, understanding the quantum behavior of electrons in semiconductor materials has made it possible to develop and produce smaller and faster semiconducting devices. Recent and promising advancements in quantum physics are now being exploited to produce new technologies. For example, certain quantum systems can be used to encode and transmit information in new ways. In addition, quantum systems can be used in optical lithographic and measuring devices to provide resolution and sensitivity limits that exceed those of non-quantum-based lithographic instruments and measuring devices.

Quantum systems comprising just two discrete states, represented by “|0

” and “|1

,” are often employed in a variety of quantum-system-based applications including quantum information encoding and processing, optical quantum lithography, and metrology. Examples of two-state quantum systems include any two photon states of an electromagnetic field, vertically and horizontally polarized photons, and the two spin states of an electron. A quantum system with two discrete states is called a “qubit system,” and the states |0

and |1

, called “qubit basis states,” can also be represented in set notation as {|0

, |1

}. A qubit system can exist in the state |0

, the state |1

, or in any of an infinite number of states that simultaneously comprise both |0

and |1

. Any of the states that include both |0

and/or |1

can be represented mathematically as a linear superposition of states: |ψ

=α|0

+β|1

The state |ψ

is called a “qubit,” and the parameters α and β are complex-valued coefficients satisfying the condition: |α|²+|β|²=1

In general, the sum of the square modulus of the coefficients is “1.” Performing a measurement on a quantum system collapses the state of the quantum system onto a basis state and produces an associated real-valued quantity. Performing a measurement on a quantum system is mathematically equivalent to projecting the state of the quantum system onto one of the basis states, and, in general, the probability of projecting the state of the quantum system onto a basis state is equal to the square of the coefficient associated with the basis state. For example, when the state |ψ

of the qubit system is measured in the basis {|0

,|1

}, one has a probability |α|² of finding the quantum system in the state |0

and a probability |β|² of finding the quantum system in the state |1

.

The infinite number of states associated with a qubit system can be geometrically represented by a unit-radius, three-dimensional sphere called a “Bloch sphere”: $\left. \psi \right\rangle = {{\cos\quad\left( \frac{\theta}{2} \right)\left. 0 \right\rangle} + {{\mathbb{e}}^{{\mathbb{i}}\quad\phi}{\sin\left( \frac{\theta}{2} \right)}\left. 1 \right\rangle}}$ where

0≦θ<π, and

0≦φ<2π.

FIGS. 1A-1C illustrate a Bloch sphere representation of a qubit system. In FIG. 1A, lines 101-103 are orthogonal x, y, and z Cartesian coordinate axes, respectively, and a Bloch sphere 106 is centered at the origin. There are an infinite number of points on the Bloch sphere 106, each point representing a unique state of a qubit system. For example, a point 108 on the Bloch sphere 106 represents a unique state of a qubit system that simultaneously comprises, in part, the state |0

and, in part, the state |1

. However, once the state of the qubit system is measured in the basis {|0

,|0

}, the state of the qubit system is projected onto the state |0

110, in FIG. 1B, or onto the state |1

112, in FIG. 1C.

The state of a combined qubit system comprising two or more qubit systems is represented by a product of qubits, each qubit associated with one of the qubit systems. For example, consider a combined qubit system comprising a first qubit system and a second qubit system that is represented by the state: |ψ

₁₂=|ψ

₁|

₂ where the state of the first qubit system is: $\left. \psi \right\rangle_{1} = {\frac{1}{\sqrt{2}}\left( {\left. 0 \right\rangle_{1} + \left. 1 \right\rangle_{1}} \right)}$ and the state of the second qubit system is: $\left. \psi \right\rangle_{2} = {\frac{1}{\sqrt{2}}\left( {\left. 0 \right\rangle_{2} + \left. 1 \right\rangle_{2}} \right)}$ The state |ψ

₂ can also be written as a linear superposition of states: $\begin{matrix} {\left. \psi \right\rangle_{12} = {\left. \psi \right\rangle_{1}\left. \psi \right\rangle_{2}}} \\ {= {\frac{1}{2}\left( {{\left. 0 \right\rangle_{1}\left. 0 \right\rangle_{2}} + {\left. 0 \right\rangle_{1}\left. 1 \right\rangle_{2}} + {\left. 1 \right\rangle_{1}\left. 0 \right\rangle_{2}} + {\left. 1 \right\rangle_{1}\left. 1 \right\rangle_{2}}} \right)}} \end{matrix}$ where the terms |0

₁, |0

₂, |0

₁|1

₂, |1

₁

₁|0

₂, and |1

₁|1

₂ are product states. Each product state in the state |ψ

₁₂ has an associated coefficient of 1/2, which indicates that when the state of the first qubit system is measured in the bases {|0

₁,|1

₁}, and the state of the second qubit system is measured in the basis {|0

₂,|1

₂}, there is a 1/4 probability of the combined qubit systems being found in any one of the product states. For example, when the states of the first and the second qubit systems are measured in the bases {|0

₁,|1

₁} and {|0

₂,|1

₂}, respectively, there is a 1/4 (|1/2|²) probability of projecting the state of the combined qubit system onto the product state |1

₁|1

₂.

The state of certain combined qubit systems cannot be represented by a product of associated qubits. These qubit systems are said to be “entangled.” Quantum entanglement is a property of quantum mechanics in which the states of two or more quantum systems are linked to one another, even though the quantum systems can be spatially separated. An example entangled state representation of an entangled two-qubit system is: $\left. \phi \right\rangle_{12} = {\frac{1}{\sqrt{2}}\left( {{\left. 0 \right\rangle_{1}\left. 1 \right\rangle_{2}} + {\left. 1 \right\rangle_{1}\left. 0 \right\rangle_{2}}} \right)}$

The entangled state |φ

₁₂ cannot be factored into a product of the qubits α₁|0

₁+β₁|1

₁ and α₂|0

₂+β₂|1

₂, for any choice of the parameters α₁, β₁, α₂, and β₂.

The state of an un-entangled, two-qubit system can be distinguished from the state of an entangled, two-qubit system as follows. Consider an un-entangled two-qubit system in the un-entangled state |ψ

₁₂. Suppose a measurement performed on the first qubit system in the basis {|0

₁,|1

₁} projects the state of the first qubit system onto the state |0

. According to the state |ψ

₁₂, the state of the un-entangled, two-qubit system just after the measurement is the linear superposition of states (|0

₁|0

₂+|0

₁|1

₂)/√{square root over (2)} When a second measurement is performed on the second qubit system in the basis {|0

₂,|1

₂} immediately following the first measurement in an identical reference frame, there is a 1/2 probability of projecting the state of the second qubit system onto the state |0

₂ and a 1/2 probability of projecting the state of the second qubit system onto the state |1

₂. In other words, the state of the second qubit system is not correlated with the state of the first qubit system. By contrast, consider an entangled, two-qubit system in the entangled state |φ

₁₂. Suppose that a first measurement performed on the first qubit system in the basis {|0

₁,|1

₁} also projects the state of the first qubit system onto the state |0

₁. According to the entangled state |φ

₁₂, the state of the entangled, two-qubit system after the first measurement is the product state |0

₁|1

₂. When a second measurement is performed on the second qubit system in the basis {|0

₂,|1

₂}, the state of the second qubit system is |1

₂ with certainty. In other words, the state of the first qubit system |0

₁ is correlated with the state of the second qubit system |1

₂.

Entangled quantum systems have a number of different and practical applications in disciplines such as metrology, cryptography, and many-qubit communications protocols. In metrology, for example, certain kinds of entangled quantum systems can be used to improve quantum-enhanced-phase measurements of a phase shift, φ, induced by an optical element on a beam of electromagnetic radiation. One conventional method for determining the phase shift φ of an optical element is to place the optical element within an internal, optical path of an interferometer. FIG. 2A illustrates an interferometer 200. The interferometer 200 comprises a first beamsplitter 202, a second beamsplitter 204, two mirrors 206 and 208, a first photon detector 210, and a second photon detector 212. The first beam splitter 202 is a 50:50 beam splitter that receives photons in input paths 214 and 216 and outputs photons in a first internal path, identified by directional arrows 218 and 220, and a second internal path, identified by directional arrows 222 and 224. The second beam splitter 204 is also a 50:50 beam splitter that recombines the photons transmitted by the first and second internal paths and outputs photons to the detectors 210 and 212. An optical element 226 with an unknown and sought after phase shift φ is located within the second internal path. The first beamsplitter 202 receives a single photon from either the input path 214 or the input path 216 and outputs a photon in the linear superposition of states: $\left. \psi \right\rangle = {\frac{1}{\sqrt{2}}\left( {\left. 0 \right\rangle + \left. 1 \right\rangle} \right)}$ where

|0

is a state that represents a photon transmitted in the first internal path, and

|1

is a state that represents a photon transmitted in the second internal path. Just before the second beamsplitter 204, the optical element 226 causes a photon in the state |1

to acquire the phase shift φ, which places the qubit in the state: $\left. \psi \right\rangle = {\frac{1}{\sqrt{2}}\left( {\left. 0 \right\rangle + {{\mathbb{e}}^{{\mathbb{i}}\quad\varphi}\left. 1 \right\rangle}} \right)}$ The phase shift φ is determined from a quantum interference pattern that is detected by the detectors 210 and 212 and is built-up over a large number of trails. The state just before the second beamsplitter 204 for N trials is a product of N qubits and is given by: $\left. \psi \right\rangle^{N} = {\prod\limits_{j = 1}^{N}{\frac{1}{2^{N/2}}\left( {\left. 0 \right\rangle_{j} + {{\mathbb{e}}^{{\mathbb{i}}\quad\varphi}\left. 1 \right\rangle_{j}}} \right)}}$ Photon interference at the second beamsplitter 204 produces over N trials a sinusoidal interference pattern of the form N cos(φ) that is detected by the detectors 210 and 212. The phase shift φ can be read from the interference pattern, and the uncertainty associated with an N trial determination of the phase shift φ, called the “shot-noise limit,” is given by: ${\Delta\quad\varphi_{SL}} = \frac{1}{\sqrt{N}}$

The expression for the shot noise limit is determined in accordance with estimation theory, as described in the book “Quantum Detection and Estimation Theory,” by C. W. Helstrom, Mathematics in Science and Engineering 123, Academic Press, New York, 1976. The shot noise limit is a value that represents the statistical dispersion, or standard deviation, in the experimentally determined value of the phase shift φ. Because the shot noise limit is inversely proportional to the square root of the number of trials, the uncertainty in the phase shift φ decreases as the number of trials N increases.

However, methods and systems that employ entangled qubit systems, such as an entangled N-qubit system in a “NOON” state, can be used to improve upon the shot noise limit by providing the optimal obtainable accuracy in the value of the phase shift φ. The N-qubit NOON state is a member of a larger class of maximally entangled qubits called “Greenberger-Horne-Zeilinger” (“GHZ”) states and is given by: $\left. \psi \right\rangle_{NOON} = {\frac{1}{\sqrt{2}}\left( {\left. {N,0} \right\rangle + \left. {0,N} \right\rangle} \right)}$ where

|N,0

represents the product state |0

₁ . . . |0

_(N), and

|0,N

represents the product state |1

₁ . . . |1

_(N). FIG. 2B illustrates use of a NOON state to determine a phase shift induced by the optical medium 226. Each state in |N,0

represents a photon transmitted in a path 230, and each state in |0,N

represents a photon transmitted in a path 232. Each photon in the state |0,N

passes through the optical medium 226 and acquires a phase shift φ, which results in the state: $\left. \psi_{N} \right\rangle = {\frac{1}{\sqrt{2}}\left( {\left. {N,0} \right\rangle + {{\mathbb{e}}^{{\mathbb{i}}\quad N\quad\varphi}\left. {0,N} \right\rangle}} \right)}$ where e^(iNφ) represents a relative phase shift difference between the photons in the state |N,0

and photons in the state |0,N

. The paths 230 and 232 intersect at a point on a substrate 234, which produces a detectable interference pattern of the form cos(Nφ) along the top surface of the substrate 234. The uncertainty associated with using the NOON state to determine the phase shift φ is, in accordance with estimation theory, given by: ${\Delta\quad\varphi_{HL}} = \frac{1}{N}$

The uncertainty Δφ_(HL) is called the “Heisenberg limit” and is the optimal accuracy obtainable for an N-trial determination of the phase shift φ. In other words, by employing N-qubit systems entangled in the NOON state, the uncertainty associated with the experimentally determined value of the phase shift φ is minimized, which is an improvement over the shot noise limit by a factor of √{square root over (N)} for the same number of trials N. However, existing methods for determining a phase shift φ using a NOON state typically employ gates and/or qubit interactions that are often cumbersome and complicated to implement. Physicists, metrologists, cryptographers, computer scientists, and quantum-information users have recognized a need for new methods that can be used to determine a phase shift φ using NOON states.

SUMMARY OF THE INVENTION

Various embodiments of the present invention are directed to methods for determining a phase shift acquired by an entangled N-qubit system represented by a NOON state. In one embodiment of the present invention, for M trials, a probe electromagnetic field may be coupled with each qubit system in the entangled N-qubit system in order to produce a probe-qubit entangled quantum system. The phase shift acquired by the qubit systems may be transferred to the probe electromagnetic field by transforming the state of each qubit-system into a linear superposition of qubit basis states. An intensity measurement may be performed on the probe electromagnetic field in order to obtain a corresponding measurement result. The measurement result falls within a measurement-result interval of a larger continuous spectrum of measurement results, and a counter associated with the measurement-result interval is incremented after each trial. For the M trials, a frequency distribution may be produced by normalizing the counter associated with each measurement-result interval. A phase shift may be determined by fitting a probability distribution associated with the probe electromagnetic field to the frequency distribution as a function of the phase shift.

BRIEF DESCRIPTION OF THE DRAWINGS

FIGS. 1A-1C illustrate a Bloch sphere representation of a qubit.

FIGS. 2A-2B illustrate two systems for determining a phase shift induced by an optical element on a beam of electromagnetic radiation.

FIG. 3 illustrates a cubic cavity.

FIG. 4 illustrates a three-dimensional coordinate system with two independent polarization vectors and a normalized wavevector as basis vectors.

FIG. 5 illustrates a representation of electric and magnetic field components of electromagnetic radiation in the coordinate system shown in FIG. 4.

FIG. 6 is an energy-level diagram of quantized electromagnetic radiation.

FIG. 7 is a phase space diagram of a number state.

FIG. 8 illustrates a probability distribution function associated with a pulse output from a source and transmitted to a detector.

FIG. 9 is a phase-space diagram of a coherent state.

FIG. 10 is a plot of a probability distribution function associated with a coherent state.

FIG. 11 illustrates a first implementation for determining a phase shift φ acquired by an entangled 4-qubit state that represents an embodiment of the present invention.

FIG. 12 is a phase space diagram of a probe field that represents an embodiment of the present invention.

FIG. 13 shows four phase-space diagrams associated with each positive phase shift θ acquired by a probe field input to the first implementation, in FIG. 11, that represent embodiments of the present invention.

FIG. 14 shows four phase-space diagrams associated with each negative phase shift θ acquired by a probe field input to the implementation, in FIG. 11, that represent embodiments of the present invention.

FIGS. 15A-15B are diagrams of example frequency distributions associated with a large number of hypothetical homodyne-detection results that represent embodiments of the present invention.

FIG. 16 is a diagram of a probability distribution fit to the frequency distribution shown in FIG. 15A that represents an embodiment of the present invention.

FIG. 17 illustrates a second implementation for determining the phase shift φ acquired by an entangled 4-qubit state that represents an embodiment of the present invention.

FIGS. 18A-18B are phase-space diagrams of a phase-shifted probe field output from a rotate quadrature operation, in FIG. 17, that represents an embodiment of the present invention.

FIG. 19 illustrates a third implementation for determining a phase shift φ acquired by an entangled 4-qubit state that represents an embodiment of the present invention.

FIG. 20 shows four phase-space diagrams associated with each increase in the average number of photons acquired by a probe field input to the third implementation, in FIG. 19, that represent an embodiment of the present invention.

FIG. 21 shows four phase-space diagrams associated with each decrease in the average number of photons acquired by a probe field input to the implementation, in FIG. 19, that represent an embodiment of the present invention.

FIG. 22 is a phase-space diagram of field states output from a displacement operation, in FIG. 19, that represents an embodiment of the present invention.

FIG. 23 is a phase-space diagram of field states output from the interaction region, in FIG. 19, that represent an embodiment of the present invention.

FIG. 24 is a phase-space diagram of field states output from the displacement operation, in FIG. 19, that represents an embodiment of the present invention.

FIG. 25 illustrates a fifth implementation for determining a phase shift φ acquired by an entangled 4-qubit state that represents an embodiment of the present invention.

FIG. 26 illustrates a sixth implementation for determining a phase shift φ acquired by an entangled 4-qubit state that represents an embodiment of the present invention.

FIG. 27 illustrates a seventh implementation comprising a single interaction region that represents an embodiment of the present invention.

DETAILED DESCRIPTION OF THE INVENTION

Various embodiments of the present invention are directed to methods for extracting a relative phase shift acquired by entangled, N-qubit systems. The N-qubit systems may be represented by a NOON state. Embodiments of the present invention are mathematical in nature and, for this reason, are described below with reference to numerous equations and numerous graphical illustrations. In order to assist in understanding descriptions of various embodiments of the present invention, an overview of quantum mechanics is provided below in a first subsection. An overview of electromagnetic radiation and quantum optics is provided in a second subsection. An overview of coherent states is provided in a third subsection. An overview of quantum entanglement is provided in a fourth subsection. Finally, various method embodiments of the present invention are described in a fifth subsection.

Overview of Quantum Mechanics

Embodiments of the present invention employ concepts in quantum mechanics. The textbook “Quantum Mechanics, Vol. I and II,” by Claude Cohen-Tannoudji, Bernard Diu and Frank Laloe, Hermann, Paris, France, 1977, is one of many references for the field of quantum mechanics. In this subsection, topics in quantum mechanics that relate to embodiments of the present invention are described. Additional details can be obtained from the above-referenced textbook, or from many other textbooks, papers, and journal articles related to quantum mechanics.

Quantum mechanics models the observed behavior, atomic and subatomic levels, of systems comprising photons, electrons, atoms, and molecules. Quantum systems exist in discrete states that are characterized by discrete measurable quantities. A state of a quantum system is represented by a ket and is denoted |Ψ

, where Ψ is a label that represents a state of a quantum system. For example, an electron has two intrinsic spin-angular-momentum states that correspond to two measurable spin-angular-momentum values

/2 and −

/2, where

is approximately 1.0546×10⁻³⁴ Js. The spin state that corresponds to the spin-angular momentum

/2 is referred to as “spin up” and is denoted |↑

, and the spin state that corresponds to the spin angular momentum −

/2 is referred to as “spin down” and is denoted |↓

. Various different labels can be assigned to various different quantum states. For example, the spin up and spin down states |↑

and |↓

can also be represented by the kets |1/2

and |−1/2

, respectively. Also, a single label can be used to represent different states in entirely different quantum systems. For example, the ket “|1

” can represent a first quantized vibrational level of a diatomic molecule and can also be used to represent a single photon, as described below, in a following subsection.

A measurement employed to determine a measurable quantity of a quantum system, such as the spin angular momentum of an electron, is represented by an operator {circumflex over (Ψ)}, where the symbol “ˆ” denotes an operator. In general, an operator operates on a ket from the left as follows: {circumflex over (Ψ)}(|Ψ

)={circumflex over (Ψ)}|Ψ

where {circumflex over (Ψ)}|Ψ

is a ket representing an observed quantum state. Typically, an operator {circumflex over (Ψ)} is associated with a set of states called “eigenstates.” An eigenstate is represented as “|ψ_(i)” with the following property: {circumflex over (Ψ)}|ψ_(i)

=ψ_(i)|ψ_(i)

where

i is a non-negative integer, and

ψ_(i) is a real value, called an “eigenvalue,” that corresponds to a discrete measurable quantity that is observed when the quantum system is in the eigenstate |ψ_(i)

.

For example, a measurement employed to determine the spin angular momentum of an electron is represented by Ŝ_(z), and the eigenvalue-eigenstate representations of observed spin-angular-momentum values are: ${{{\hat{S}}_{z}\left. \uparrow \right\rangle} = {\frac{\hslash}{2}\left. \uparrow \right\rangle}},{and}$ ${{\hat{S}}_{z}\left. \downarrow \right\rangle} = {{- \frac{\hslash}{2}}\left. \downarrow \right\rangle}$

The eigenstates are basis vectors of a complex vector space called a “Hilbert space,” and the number of eigenstates is the dimension of the Hilbert space. For example, a Hilbert space of an electron is two-dimensional, with eigenstates |↑

and |↓

. A Hilbert space with N eigenstates {|ψ_(i)

} is N-dimensional, and any state |Ψ

in the Hilbert space can be written as a linear superposition of the eigenstates as follows: $\left. \Psi \right\rangle = {\sum\limits_{i = 1}^{N}{c_{i}\left. \psi_{i} \right\rangle}}$ where c_(i) is a complex valued coefficient called the “amplitude.” A Hilbert space also includes a mathematical operation called the “inner product.” The inner product of two states |Ψ

and |Ξ

is represented by:

Ξ|Ψ

where

Ξ| is called a “bra,” and represents the complex conjugate and transpose of the state |Ξ

. The inner product has the following property:

Ξ|Ψ

=

Ψ|Ξ

* where “*” represents the complex conjugate. The basis eigenstates of a Hilbert space are orthonormal, or in mathematical notation:

ψ_(i)|ψ_(j)

=δ_(ij) where δ_(ij) is 1 when i equals j, and 0 otherwise. For example, the inner product of the eigenstates of a single electron Hilbert space are:

↑|↑

=

↓|↓

=1, and

↑|↓

=

↓|↑

=0

The orthonomality property of the eigenstates of a Hilbert space can be used to determine the coefficients of the linear superposition of states |Ψ

. Taking the inner product of |Ψ

with

ψ_(j)| gives the corresponding coefficient: $\left\langle \psi_{j} \middle| \Psi \right\rangle = {{\sum\limits_{i = 1}^{N}{c_{i}\left\langle \psi_{j} \middle| \psi_{i} \right\rangle}} = {{\sum\limits_{i = 1}^{N}{c_{i}\delta_{ij}}} = c_{j}}}$ Substituting for the coefficients in the linear superposition gives: $\left. \Psi \right\rangle = {\sum\limits_{i = 1}^{N}{\left. \psi_{i} \right\rangle\left\langle \psi_{i} \middle| \Psi \right\rangle}}$ Because |Ψ

is an arbitrary ket in the Hilbert space, ${\sum\limits_{i = 1}^{N}{\left. \psi_{i} \right\rangle\left\langle \psi_{i} \right.}} = \hat{1}$ where “Î” is the identity operator. The summation is called the “completeness relation,” and the eigenstates {|ψ_(i)} are said to be “complete.”

Because eigenstates of a Hilbert space are orthonormal and provide a basis for the Hilbert space, the eigenstates can be represented by orthogonal normalized column vectors and an operator can be represented by square matrix. For example, the eigenstates of a single electron Hilbert space are represented by the column vectors: ${\left. \uparrow \right\rangle \doteq \begin{bmatrix} 1 \\ 0 \end{bmatrix}},{{{and}\quad\left. \downarrow \right\rangle} \doteq \begin{bmatrix} 0 \\ 1 \end{bmatrix}}$ where the symbol “{dot over (=)}” stands for “is represented by.” The transposed complex conjugates of the eigenstates are represented by the row vectors:

↑|{dot over (=)}[1 0],and

↓|{dot over (=)}[0 1] Using the completeness relation, an operator Ô on the basis {|ψ_(i)

} can also be represented by: $\hat{O} = {\sum\limits_{i = 1}^{N}{\sum\limits_{j = 1}^{N}{\left. \psi_{i} \right\rangle\left\langle {\psi_{i}{\hat{O}}\psi_{j}} \right\rangle\left\langle \psi_{j} \right.}}}$ where

ψ_(i)|Ô|ψ_(j)

is a matrix element. The matrix corresponding to the operator Ô on the basis {|ψ_(i)

} can be represented as follows: $\hat{O} \doteq \begin{bmatrix} \left\langle {\psi_{1}{\hat{O}}\psi_{1}} \right\rangle & \left\langle {\psi_{1}{\hat{O}}\psi_{2}} \right\rangle & \ldots & \left\langle {\psi_{1}{\hat{O}}\psi_{N}} \right\rangle \\ \left\langle {\psi_{2}{\hat{O}}\psi_{1}} \right\rangle & \left\langle {\psi_{2}{\hat{O}}\psi_{2}} \right\rangle & \quad & \vdots \\ \vdots & \quad & ⋰ & \quad \\ \left\langle {\psi_{N}{\hat{O}}\psi_{1}} \right\rangle & \ldots & \quad & \left\langle {\psi_{N}{\hat{O}}\psi_{N}} \right\rangle \end{bmatrix}$ The matrix representation of the operator Ô equal to {circumflex over (Ψ)} has zero off diagonal elements, and the diagonal elements are the eigenvalues {ψ_(i)}. For example, the electron spin operator can be given by: ${\hat{S}}_{z} = {\frac{\hslash}{2}{\hat{\sigma}}_{z}}$ where {circumflex over (σ)}_(z)=|↑

↑|−|↓

↓|. The matrix representation of the electron spin operator Ŝ_(z) is given by: $\begin{matrix} {{\hat{S}}_{z}\overset{.}{=}\begin{bmatrix} \left\langle \uparrow{{\hat{S}}_{z}}\uparrow \right\rangle & \left\langle \uparrow{{\hat{S}}_{z}}\downarrow \right\rangle \\ \left\langle \downarrow{{\hat{S}}_{z}}\uparrow \right\rangle & \left\langle \downarrow{{\hat{S}}_{z}}\downarrow \right\rangle \end{bmatrix}} \\ {= {\frac{\hslash}{2}\begin{bmatrix} 1 & 0 \\ 0 & {- 1} \end{bmatrix}}} \end{matrix}$ An operator Ô is called a “Hermitian operator” if Ô=Ô^(†) The corresponding matrix elements satisfies the condition:

ψ_(i)|Ô|ψ_(j)

=

ψ_(j)|Ô^(†)|ψ_(i)

Prior to a measurement, a quantum system can simultaneously exist in all of the eigenstates of a corresponding Hilbert space, which is represented by the (pure state) linear superposition of states: $\left. \Psi \right\rangle = {{\sum\limits_{i = 1}^{N}{c_{i}\left. \psi_{i} \right\rangle}} = {\sum\limits_{i = 1}^{N}{\left. \psi_{i} \right\rangle\left\langle \psi_{i} \middle| \Psi \right\rangle}}}$ A measurement performed on the quantum system projects the state |Ψ

of the quantum system onto one of the eigenstates. In other words, a measurement on a quantum system is essentially a filtering process that places the quantum system in one of the eigenstates in the linear superposition at the time of the measurement. For example, an electron with an unknown spin orientation prior to a measurement is represented in a linear superposition of states: |Ψ

=c ₁|↑

+c ₂↓

A spin determination measurement Ŝ_(z) projects the state of the electron onto either the state |↑

or the state |↓

at the time of the measurement. In other words, just after the spin determination measurement, the electron is either in the state |↑

or the state |↓

.

There is a corresponding irreversible change to the state of a quantum system as a result of a measurement. Irreversibility can only be avoided when the quantum system is already in one of the quantum states before the measurement is performed. As a result, one cannot infer the prior state of a quantum system based on the outcome of a single measurement. For example, if the outcome of a spin measurement is

/2, it is not possible to determine whether the system was already in the state |↑

or in a linear superposition of the spin states |↑

and |↓

at the time of the measurement.

Although it is not possible to know in advance which of the various states |ψ_(i)

the state of a quantum system will be projected onto, the probability of the quantum system being found in a particular state |ψ_(i)

immediately after the measurement is given by:

Probability for ψ_(i)=|c_(i)|²=|

ψ_(i)|Ψ

|² where |Ψ

is normalized, and |c_(i)|² equals c_(i)*c_(i) and gives the outcome probability. For example, prior to a spin determination measurement in the spin basis {|↑

,|↓

}, consider an electron coherently prepared with a 1/2 probability of being found in the spin state |↑

and a 1/2 probability of being found in the spin state |↓

. The state associated with the electron in such as spine state prior to a spin determination measurement can be represented by: $\left. \Psi \right\rangle = {{\frac{1}{\sqrt{2}}\left. \uparrow \right\rangle} + {\frac{1}{\sqrt{2}}\left. \downarrow \right\rangle}}$

The expectation value of a measurement performed on an ensemble of quantum systems that are represented by the linear superposition of states |Ψ

is mathematically represented by:

{circumflex over (Ψ)}

=

Ψ|{circumflex over (Ψ)}|Ψ

and is determined by applying the completeness relation as follows: $\begin{matrix} {\left\langle \hat{\Psi} \right\rangle = {\sum\limits_{i = 1}^{N}{\sum\limits_{j = 1}^{N}{\left\langle \Psi \middle| \psi_{i} \right\rangle\left\langle {\psi_{i}{\hat{\Psi}}\psi_{j}} \right\rangle\left\langle \psi_{j} \middle| \Psi \right\rangle}}}} \\ {= {\sum\limits_{i = 1}^{N}{\psi_{i}{\left\langle \psi_{i} \middle| \Psi \right\rangle }^{2}}}} \end{matrix}$ The expectation value represents the weighted eigenvalue average result expected from measurements on the quantum systems in the ensemble, where the initial state |Ψ

of the quantum system is the same for each member of the ensemble. In other words, the linear superposition of states representing each quantum system is identical prior to the measurement. In practice, such an ensemble could be realized by preparing many identical and independent quantum systems all in the same state, or by repeatedly preparing a single system in the same state. Note that the expectation value may not be the value obtained for each measurement and, therefore, is not to be confused with the eigenvalue obtained from the measurement. For example, the expectation value of Ŝ₂ can be any real value between the eigenvalues

/2 and

/2, but the actual measured value of S, for an electron is always either

/2 or −

/2 in each individual measurement.

A tensor product is a way of combining Hilbert spaces of different quantum systems to form Hilbert spaces that represent combined quantum systems. For example, H_(Ψ) is a Hilbert space of a first quantum system, and H_(Ξ) is a Hilbert space of a second quantum system. The Hilbert space denoted by H_(Ψ)

H_(Ξ) represents a combined Hilbert space, where the symbol

represents a tensor product. The operators {circumflex over (Ψ)} and {dot over (Ξ)} correspond to the Hilbert spaces H_(Ψ) and H_(Ξ), respectively, and each operates only on the corresponding eigenstates as follows: ({circumflex over (ψ)}

{circumflex over (Ξ)})(|ψ

ξ

)=(Ψ|ψ

)

({circumflex over (Ξ)}|ξ

) where |ψ

represents a state in the Hilbert space H_(Ψ), and |ξ

represents a state in the Hilbert space H_(Ξ). The tensor product |ψ

ξ

can be abbreviated as |ψ

|ξ

, |ψ,ξ

or |ψξ

. The spin states of two electrons in an atomic orbital is an example of a combined Hilbert space. The two electrons can either both be spin up, both be spin down, the first electron spin up and the second electron spin down, or the first electron spin down and the second electron spin up. The various tensor product representations of two spin up electrons are given by: |↑

₁

|↑

₂=|↑

₁|↑

₂=|↑,↑

₁₂=|↑↑

₁₂ where the subscripts 1 and 2 refer to the first and second electrons.

In quantum mechanics, there are also measurable quantities with continuous eigenvalue spectrums. The dimensionality of the corresponding Hilbert spaces are infinite and many of the properties described above for discrete quantum systems can be generalized for continuous quantum systems. A continuous eigenvalue equation is: {circumflex over (ζ)}|ζ

=ζ|ζ

where ζ represents a continuous eigenvalue, and the ket |ζ

is a continuous eigenstate of the operator {circumflex over (ζ)}. For example, for an unbound particle in one dimension, both position q and momentum p are continuous eigenvalues of the position and momentum operators {circumflex over (q)} and {circumflex over (p)}, respectively, and can assume any real value between −∞ and ∞.

The properties of the continuous variable ζ can be generalized as follows: ⟨ζ|ζ^(′)⟩ = δ(ζ − ζ^(′)), ∫_(−∞)^(∞)𝕕ζζ⟩⟨ζ = 1, and ⟨ζζ̂ζ^(′)⟩ = ζ^(′)δ(ζ − ζ^(′)), where δ(ζ−ζ′) is the delta function, which has numerous limit representations, such as ${\delta\left( {\zeta - \zeta^{\prime}} \right)} = {\lim\limits_{\Delta->0}{\frac{1}{\sqrt{2\quad\pi\quad\Delta^{2}}}{\exp\left( {- \frac{\left( {\zeta - \zeta^{\prime}} \right)^{2}}{2\quad\Delta^{2}}} \right)}}}$ A state ket for an arbitrary physical state can be expanded in terms of the states {|ζ

} as follows: α⟩ = ∫_(−∞)^(∞)𝕕ζζ⟩⟨ζ|α⟩ For example, consider placing in the path of a particle a detector that outputs the position of the particle when the particle is at the position q. Immediately after the measurement is taken, the system, initially in the state |α

, is projected onto the state represented by |q

in much the same way an arbitrary electron-spin state is projected onto one of the two spin states when a spin detection measurement is performed. Other properties of the continuous variable ζ are given by: ∫_(−∞)^(∞)𝕕ζ⟨ζ|α⟩² = 1, and ⟨β|α⟩ = ∫_(−∞)^(∞)𝕕ζ⟨β|ζ⟩⟨ζ|α⟩

The momentum operator {circumflex over (p)} can also be represented by a differential operator −i

∂/∂q. As a result, both the position and momentum operators satisfy the canonical commutation relations: [{circumflex over (q)}_(i),{circumflex over (q)}_(j)]=0 [{circumflex over (p)}_(i),{circumflex over (p)}_(j)]=0, and [{circumflex over (q)}_(i),{circumflex over (p)}_(j)]=ihδ_(ij) where

i and j represent orthogonal coordinates, such as the Cartesian x, y, and z coordinates, and

the commutator is defined as [A, B]=AB−BA.

An Overview of Electromagnetic Radiation and Quantum Optics

In this subsection, a brief description of electromagnetic radiation and quantum optics that relates to embodiments of the present invention is described. The textbooks “Quantum Optics,” M. O. Scully and M. S. Zubairy, Cambridge University Press, Cambridge, United Kingdom, 1997, and “The Quantum Theory of Light (3^(rd) Edition),” R. Loudon, Oxford University Press, New York, 2000 are two of many references for quantum optics. Additional details can be obtained from the above-referenced textbooks, or from many other textbooks, papers, and journal articles in this field.

Quantum optics is a field of physics that relates the application of quantum mechanics to electromagnetic radiation. Electromagnetic radiation confined to a cavity with perfectly reflecting walls is quantized. Quantized electromagnetic radiation can be applied to more general unconfined optical systems, such as electromagnetic radiation propagating in free space or in an optical fiber.

Electromagnetic radiation confined to a cavity, with no free charges and currents, comprises an electric field component {right arrow over (E)}({right arrow over (r)},t) and a magnetic field component {right arrow over (B)}({right arrow over (r)},t) that are related in terms of a vector potential {right arrow over (A)}({right arrow over (r)},t) satisfying the wave equation: ${{\nabla^{2}\overset{\_}{A}} - {\frac{1}{c^{2}}\frac{\partial^{2}\overset{\_}{A}}{\partial t^{2}}}} = 0$ and the Coulomb, non-relativistic gauge condition: ∇·{right arrow over (A)}({right arrow over (r)},t)=0 where the electric and magnetic field components are determined by: ${{\overset{\_}{E}\left( {\overset{\_}{r},t} \right)} = {- \frac{\partial{\overset{\_}{A}\left( {\overset{\_}{r},t} \right)}}{\partial t}}},{and}$ ${\overset{\_}{B}\left( {\overset{\_}{r},t} \right)} = {\nabla{\times {\overset{\_}{A}\left( {\overset{\_}{r},t} \right)}}}$ The electromagnetic radiation is assumed to be propagating but is subject to periodic boundary conditions imposed by a cubic, or quantization, cavity with perfectly reflecting walls, where the walls are of length L. FIG. 3 illustrates a cubic cavity 300. Orthogonal axes 302, 304, and 306 represent the x, y, and z Cartesian coordinate axes. The finite dimensional cubic cavity 300 imposes periodic boundary conditions on solutions to the wave equation. For example, in the x, y, and z-directions, plane wave solutions to the vector potential wave equation satisfy the condition: exp(i{right arrow over (k)}·{right arrow over (r)})=exp( i{right arrow over (k)}·({right arrow over (r)}+ L )) where L is vector (L, L, L), and

{right arrow over (k)} is called the “wavevector” with components: ${\overset{\_}{k} = {\frac{2\quad\pi}{L}\left( {m_{x},m_{y},m_{z}} \right)}},{and}$

m_(x), m_(y), and m_(z) are integers.

Each set of integers (m_(x), m_(y), m_(z)) specifies a normal mode of the electromagnetic radiation, and the magnitude of the wavevector {right arrow over (k)}, k, is equal to ω_(k)/c, where c represents the speed of light in free space and ω_(k) is the angular frequency. Note that in real life the spectrum of normal modes of an electromagnetic field is actually continuous and a discrete spectrum of normal modes suggested by the wavevector {right arrow over (k)} is an approximation to the continuous spectrum.

A propagating vector potential solution to the wave equation above that satisfies the periodic boundary conditions is: ${A\left( {r,t} \right)} = {\sum\limits_{\overset{\_}{k},s}{{\overset{\_}{e}}_{\overset{\_}{k}s}\left( {{A_{\overset{\_}{k}s}{\mathbb{e}}^{{\mathbb{i}}\quad{({{\overset{\_}{k} \cdot \overset{\_}{r}} - {\omega_{k}t}})}}} + {A_{\overset{\_}{k}s}^{*}{\mathbb{e}}^{{- {\mathbb{i}}}\quad{({{\overset{\_}{k} \cdot \overset{\_}{r}} - {\omega_{k}t}})}}}} \right)}}$ where

A_({right arrow over (k)}s) is a complex amplitude of the electromagnetic radiation;

{right arrow over (e)}_({right arrow over (k)}s) represents two unit-length polarization vectors; and

m_(x), m_(y), m_(z)=0, ±1, ±2, ±3 . . .

The sum over {right arrow over (k)} represents the sum over the integers (m_(x), m_(y), m_(z)), and the sum over s is the sum over the two independent polarizations that are associated with each {right arrow over (k)}. The two polarization vectors are orthogonal as indicated by: {right arrow over (e)} _({right arrow over (k)}s) ·{right arrow over (e)} _({right arrow over (k)}s′)=δ_(ss′), and from the gauge condition given above: {right arrow over (k)}·{right arrow over (e)} _({right arrow over (k)}s)=0, for both polarization directions s. The two polarization vectors {right arrow over (e)}_({right arrow over (k)}1) and {right arrow over (e)}_({right arrow over (k)}2) form a right-handed coordinate system with a normalized wavevector given by: ${{\overset{\_}{e}}_{\overset{\_}{k}1} \times {\overset{\_}{e}}_{\overset{\_}{k}2}} = {\frac{\overset{\_}{k}}{\overset{\_}{k}} = \overset{\_}{\kappa}}$ FIG. 4 illustrates a three-dimensional right-handed coordinate system with two independent polarization vectors {right arrow over (e)}_({right arrow over (k)}s) and a normalized wavevector {right arrow over (κ)} as basis vectors. In FIG. 4, the wavevector {right arrow over (k)} 402, and the polarization vectors, {right arrow over (e)}_({right arrow over (k)}1) 404 and {right arrow over (e)}_({right arrow over (k)}2) 406, are three orthogonal unit length basis vectors of a coordinate system with coordinate axes represented by lines 408, 410, and 412, respectively.

The propagating electric and magnetic field components of the vector potential are: ${{\overset{\_}{E}\left( {\overset{\_}{r},t} \right)} = {{\mathbb{i}}\quad{\sum\limits_{\overset{\_}{k},s}{\omega_{k}{{\overset{\_}{e}}_{\overset{\_}{k}s}\left\lbrack {{A_{\overset{\_}{k}s}{\mathbb{e}}^{{\mathbb{i}}\quad{({{\overset{\_}{k} \cdot \overset{\_}{r}} - {\omega_{k}t}})}}} - {A_{\overset{\_}{k}s}^{*}{\mathbb{e}}^{{- {\mathbb{i}}}\quad{({{\overset{\_}{k} \cdot \overset{\_}{r}} - {\omega_{k}t}})}}}} \right\rbrack}}}}},{and}$ ${\overset{\_}{B}\left( {\overset{\_}{r},t} \right)} = {\frac{\mathbb{i}}{c}{\sum\limits_{\overset{\_}{k},s}{{\omega_{k}\left( {\overset{\_}{\kappa} \times {\overset{\_}{e}}_{\overset{\_}{k}s}} \right)}\left\lbrack {{A_{\overset{\_}{k}s}{\mathbb{e}}^{{\mathbb{i}}\quad{({{\overset{\_}{k} \cdot \overset{\_}{r}} - {\omega_{k}t}})}}} - {A_{\overset{\_}{k}s}^{*}{\mathbb{e}}^{{- {\mathbb{i}}}\quad{({{\overset{\_}{k} \cdot \overset{\_}{r}} - {\omega_{k}t}})}}}} \right\rbrack}}}$ Both the electric field {right arrow over (E)}({right arrow over (F)},t) and magnetic field {right arrow over (B)}({right arrow over (r)},t) are propagating wave solutions referred to as the “classical” representation of the electric and magnetic field, are orthogonal to one another, and are both orthogonal to the wavevector {right arrow over (k)}.

FIG. 5 illustrates a representation of electric and magnetic field components of electromagnetic radiation in the right-handed coordinate system shown in FIG. 4. The electromagnetic radiation is directed along the wavevector {right arrow over (κ)} 402 axis. The electric field component {right arrow over (E)}({right arrow over (r)},t) 502 and magnetic field component {right arrow over (B)}({right arrow over (r)},t) 504 are directed along the orthogonal polarization vectors {right arrow over (e)}_({right arrow over (k)}1) 404 and {right arrow over (e)}_({right arrow over (k)}2) 406, respectively, and appear frozen at a particular time t.

The energy of the electromagnetic radiation can be determined by evaluating the Hamiltonian: $\begin{matrix} {H = {\frac{1}{2}{\int_{V}{\left( {{ɛ_{0}{\overset{\_}{E} \cdot \overset{\_}{E}}} + {\frac{1}{\mu_{0}}{\overset{\_}{B} \cdot \overset{\_}{B}}}} \right){\mathbb{d}V}}}}} \\ {{= {2\quad ɛ_{0}V{\sum\limits_{\overset{\_}{k},s}{\omega_{k}^{2}A_{\overset{\_}{k}s}A_{\overset{\_}{k}s}^{*}}}}},} \end{matrix}$ where

ε₀ is the electric permittivity of free space,

μ₀ is the magnetic permeability of free space, and

V is the volume of the cavity.

The electric permittivity ε₀ represents the degree to which a vacuum space can store electrical potential energy under the influence of an electric field, and the magnetic permeability μ₀ represents the degree to which the vacuum modifies the flux of a magnetic field. In a non-conducting medium, the electric permittivity is further multiplied by ε, which is the degree to which the medium enhances the storage of electrical potential energy, and the magnetic permeability is further multiplied by μ, which is the degree to which the medium further enhances the flux of a magnetic field.

In order to quantize the electric field {right arrow over (E)}({right arrow over (r)},t) and magnetic field {right arrow over (B)}({right arrow over (r)},t) components, the canonical variables for position, q_({right arrow over (k)}s), and momentum, p_({right arrow over (k)}s), are introduced into the Hamiltonian by setting: $A_{\overset{\_}{k}s} = {\frac{1}{2\quad\omega_{k}\sqrt{ɛ_{0}V}}\left( {{\omega_{k}q_{\overset{\_}{k}s}} + {{\mathbb{i}}\quad p_{\overset{\_}{k}s}}} \right)}$ As a result, the Hamiltonian for the electromagnetic radiation becomes: $H = {\frac{1}{2}{\sum\limits_{\overset{\_}{k},s}\left( {p_{\overset{\_}{k}s}^{2} + {\omega_{k}^{2}q_{\overset{\_}{k}s}^{2}}} \right)}}$ Each term in the Hamiltonian is the energy of a harmonic oscillator with vibrational mode {right arrow over (k)}s, where the term p_({right arrow over (k)}s) ²/2 is the kinetic energy, and the term ω_(k)q_({right arrow over (k)}s) ²/2 is the potential energy of a harmonic oscillator with a unit mass. The Hamiltonian is quantized by replacing the position and momentum variables q_({right arrow over (k)}s) and p_({right arrow over (k)}s) with quantum mechanical position and momentum operators {circumflex over (q)}_({right arrow over (k)}s) and {circumflex over (p)}_({right arrow over (k)}s), respectively, to give the quantum Hamiltonian operator: $\hat{H} = {\frac{1}{2}{\sum\limits_{\overset{\_}{k},s}\left( {{\hat{p}}_{\overset{\_}{k}s}^{2} + {\omega_{k}^{2}{\hat{q}}_{\overset{\_}{k}s}^{2}}} \right)}}$ Annihilation and creation operators are defined by: ${{\hat{a}}_{\overset{\_}{k}s} = {\frac{1}{\sqrt{2\quad\hslash\quad\omega_{k}}}\left( {{\omega_{k}{\hat{q}}_{\overset{\_}{k}s}} + {{\mathbb{i}}\quad{\hat{p}}_{\overset{\_}{k}s}}} \right)}},{and}$ ${{\hat{a}}_{\overset{\_}{k}s}^{\dagger} = {\frac{1}{\sqrt{2\quad\hslash\quad\omega_{k}}}\left( {{\omega_{k}{\hat{q}}_{\overset{\_}{k}s}} - {{\mathbb{i}}\quad{\hat{p}}_{\overset{\_}{k}s}}} \right)}},$ and substituting the annihilation and creation operators in the quantum Hamiltonian operator gives: ${\hat{H} = {\sum\limits_{\overset{\_}{k},s}{{\hslash\omega}_{k}\left( {{{\hat{a}}_{\overset{\_}{k}s}^{\dagger}{\hat{a}}_{\overset{\_}{k}s}} + \frac{1}{2}} \right)}}},$ where {right arrow over (a)}_({right arrow over (k)}s) ^(†)â_({right arrow over (k)}s) is called the “number operator” and is also denoted by {circumflex over (n)}_({right arrow over (k)}s). Using the canonical commutation relations for the position and momentum operators, the annihilation and creation operators satisfy the commutation relations given by: [â_({right arrow over (k)}s),â_({right arrow over (k)}′s′]=)0=[â_({right arrow over (k)}s) ^(†),â_({right arrow over (k)}′s′) ^(†)], and [â_({right arrow over (k)}s),â_({right arrow over (k)}′s′) ^(†)]=δ_({right arrow over (kk′)})δ_(ss′)

When the electromagnetic radiation is quantized, the amplitudes A_({right arrow over (k)}s) become operators: ${{\hat{A}}_{\overset{\_}{k}s} = {\sqrt{\frac{\hslash}{2\omega_{k}ɛ_{0}V}}{\hat{a}}_{\overset{\_}{k}s}}},$ which can be substituted in the classical electric and magnetic field equations above to obtain electric and magnetic field operators: ${{\hat{E}\left( {\overset{\_}{r},t} \right)} = {{\mathbb{i}}{\sum\limits_{\overset{\_}{k},s}{\sqrt{\frac{\hslash\omega}{2ɛ_{0}V}}{e_{\overset{\_}{k}s}\left\lbrack {{{\hat{a}}_{\overset{\_}{k}s}{\mathbb{e}}^{{\mathbb{i}}{({{k \cdot r} - {\omega_{k}t}})}}} - {{\hat{a}}_{\overset{\_}{k}s}^{\dagger}{\mathbb{e}}^{- {{\mathbb{i}}{({{k \cdot r} - {\omega_{k}t}})}}}}} \right\rbrack}}}}},{and}$ ${\hat{B}\left( {\overset{\_}{r},t} \right)} = {\frac{\mathbb{i}}{c}{\sum\limits_{\overset{\_}{k},s}{\left( {\overset{\_}{\kappa} \times e_{\overset{\_}{k}s}} \right)\sqrt{\frac{\hslash\omega}{2ɛ_{0}V}}{e_{\overset{\_}{k}s}\left\lbrack {{{\hat{a}}_{\overset{\_}{k}s}{\mathbb{e}}^{{\mathbb{i}}{({{k \cdot r} - {\omega_{k}t}})}}} - {{\hat{a}}_{\overset{\_}{k}s}^{\dagger}{\mathbb{e}}^{- {{\mathbb{i}}{({{k \cdot r} - {\omega_{k}t}})}}}}} \right\rbrack}}}}$ Both the electric and magnetic field operators are Hermitian and represent measurable electric and magnetic fields.

The electric field accounts for most of the interactions with charged matter, because the magnitude of the magnetic field is smaller than the electric field by the factor 1/c. As a result, the electric field alone is generally used to characterize the behavior of electromagnetic radiation and any interactions with charged matter, and the magnetic field component can be ignored.

Quantum computation and quantum information processing systems can be operated using a single-mode {right arrow over (k)}s of electromagnetic radiation. As a result, the Hamiltonian operator for a single-mode of electromagnetic radiation reduces to: ${\hat{H} = {{\hslash\omega}\left( {{{\hat{a}}^{\dagger}\hat{a}} + \frac{1}{2}} \right)}},$ where â and â^(†) replace the mode-dependent operators â_({right arrow over (k)}) _(j) _(s) _(j) and â_({right arrow over (k)}) _(j) _(s) _(j) ^(†) in the Hamiltonian above. The eigenstates and the corresponding energy eigenvalues of the single-mode Hamiltonian are: ${{\hat{H}\left. n \right\rangle} = {{{{\hslash\omega}\left( {{{\hat{a}}^{\dagger}\hat{a}} + \frac{1}{2}} \right)}\left. n \right\rangle} = {E_{n}\left. n \right\rangle}}},$ where |n

is called a “number state,” n is a nonnegative integer called a “photon number,” and E_(n) is an energy eigenvalue.

The annihilation and creation operators operate on a number state as follows: â|n

=√{right arrow over (n)}|n−1

, â ^(†) |n

=√{square root over (n+1)}|n+1

, and ĥ|n

=n|n

, where {circumflex over (n)} represents the operator â^(†)â and is called the “number operator.” The number states can be generated by repeated application of the annihilation and creation operators to the number states. For example, repeated application of the annihilation operator to a number state lowers the photon number: ${\left. 0 \right\rangle = {\frac{{\hat{a}}^{n}}{\sqrt{n!}}\left. n \right\rangle}},$ where |0

is called the “vacuum state” and represents the lowest energy state of the electromagnetic radiation. Beginning with the vacuum state, and repeatedly applying the creation operator gives: $\left. n \right\rangle = {\frac{\left( {\hat{a}}^{\dagger} \right)^{n}}{\sqrt{n!}}\left. 0 \right\rangle}$ The number states are orthogonal and form a compete set represented by: ⟨n^(′)|n⟩ = δ_(n^(′)n), and ${\sum\limits_{n = 0}^{\infty}{\left. n \right\rangle\left\langle n \right.}} = 1$

In general, the energy eigenvalue equation associated with a number state |n

is: ${\hat{H}\left. n \right\rangle} = {{{{\hslash\omega}\left( {n + \frac{1}{2}} \right)}\left. n \right\rangle} = {E_{n}{\left. n \right\rangle.}}}$ Applying the annihilation and creation operators to the energy eigenvalue equation gives: ${{\hat{H}\left( {\hat{a}\left. n \right\rangle} \right)} = {{{{\hslash\omega}\left( {n - \frac{1}{2}} \right)}\left. {n - 1} \right\rangle} = {\left( {E_{n} - {\hslash\omega}} \right)\left. {n - 1} \right\rangle}}},{and}$ ${{\hat{H}\left( {{\hat{a}}^{\dagger}\left. n \right\rangle} \right)} = {{{{\hslash\omega}\left( {n + \frac{3}{2}} \right)}\left. {n + 1} \right\rangle} = {\left( {E_{n} + {\hslash\omega}} \right)\left. {n + 1} \right\rangle}}},$ which shows that the energy levels of electromagnetic radiation are equally spaced by a quantum of energy

ω. In other words, the excitations of electromagnetic radiation occur in discrete amounts of energy

ω called “photons.” The photon number n refers to the number of photons has comprising the electromagnetic radiation.

FIG. 6 is an energy level diagram of quantized electromagnetic radiation. Horizontal lines, such as horizontal line 602, represent energy levels of electromagnetic radiation. Energy level 604 is the lowest energy level, which corresponds to the vacuum state |0

. The energy of the vacuum state is

ω/2 or 1/2 the energy of a single photon. Higher energy levels of electromagnetic radiation are each separated by the same quantum of energy hω. For example, the energy level 606 represents electromagnetic radiation with a total electromagnetic energy of 5

ω/2, which can be thought of as the energy of two photons plus the vacuum state energy

ω/2. The annihilation operator corresponds to removal of a photon from the electromagnetic radiation, and the creation operator corresponds to addition of a photon to the electromagnetic radiation. For example, the annihilation operator â represents an electromagnetic-radiation transition 610 from the state |n

602 to the lower energy state |n−1

608. The transition 610 is achieved by giving up a photon to the surroundings. By contrast, the creation operator â^(†) represents an electromagnetic-radiation transition 614 from the state |n

602 to the higher energy state |n+1

612. The transition 614 is achieved by accepting a photon from the surroundings. Note that typically the surroundings can be an atom, a quantum dot, or any other system that couples to the field through a dipole interaction. Loss or absorption of a photon will involve a simultaneous excitation of the surrounding system and creation or emission of a photon will involve a corresponding de-excitation of the surrounding system.

Both the creation and annihilation operators are not Hermitian. As a result, the operators â and â^(†) cannot represent measurable quantities. However, the annihilation and creation operators can be used to construct the following Hermitian quadrature operators: ${\hat{X} = {{\sqrt{\frac{m\quad\omega}{2\hslash}}\hat{q}} = \left( {\hat{a} + {\hat{a}}^{\dagger}} \right)}},{and}$ $\hat{Y} = {\frac{\hat{p}}{\sqrt{2m\quad{\hslash\omega}}} = {{\mathbb{i}}\left( {{\hat{a}}^{\dagger} - \hat{a}} \right)}}$ The quadrature operators are essentially dimensionless position and momentum operators and are associated with the electric field amplitudes oscillating out of phase with each other by 90°. The energy eigenvalue can be rewritten in terms of the quadrature operators as: $\hat{H} = {{{{\hslash\omega}\left( {{\hat{X}}^{2} + {\hat{Y}}^{2}} \right)}\left. n \right\rangle} = {{{\hslash\omega}\left( {n + \frac{1}{2}} \right)}{\left. n \right\rangle.}}}$ The number states have the quadrature-operator eigenvalue property: ${{\left( {{\hat{X}}^{2} + {\hat{Y}}^{2}} \right)\left. n \right\rangle} = {\left( {n + \frac{1}{2}} \right)\left. n \right\rangle}},$ and the number states have identical properties for the {circumflex over (X)} and Ŷ quadrature operators. For example, the quadrature-operator-expectation values:

n|{circumflex over (X)}|n

=

n|Ŷ|n

=0 The quadrature operators can be used to construct a phase-space diagram of the number states.

FIG. 7 is a phase space diagram of a number state. Horizontal line 702 is the X quadrature axis, and vertical line 704 is the Y quadrature axis. The origin 706 corresponds to the vacuum state |0

, and circle 708 is a qualitative representation of a number state |n

, where the radius of the circle 708 is the number of photons n. The number state |n

can be located anywhere on the circle 708, which indicates that uncertainty in the number of photons is zero, but the phase angle associated with the number state is entirely undefined, so measurement will give all possible outcomes with equal probability.

The number states also have the property: =Σ=

n|Ê(t)|n

=0, where Σ is called a “coherent signal.” The zero valued coherent signal of a photon state is consistent with the sinusoidal variation of the electric field with time at a fixed observation point.

Photons can be generated by a photon source and transmitted through free space or in an optical fiber. The photon source can be a pulsed laser that generates a single pulse or a train of pulses, each pulse containing one or more photons that all have the same optical properties, such as wavelength and direction. Photons with the same optical properties are called “coherent.” However, the source, the detector, and a medium, such as an optical fiber, separating the source from the detector do not define an optical cavity. The source and the detector are parts of a continuous unidirectional flow of optical energy with no significant reflection or recycling of the optical energy. A pulse transmitted through free space or an optical fiber is described by a wavepacket that can be represented by a time-dependent, Gaussian-shaped function given by: ${{\xi(t)} = {\left( \frac{2\Delta^{2}}{\pi} \right)^{1/4}\exp\left\{ {{{- {\mathbb{i}\omega}_{0}}t} - {\Delta^{2}\left( {t_{0} - t} \right)}^{2}} \right\}}},$ where

ω₀ is the central frequency of the pulse spectrum,

t is time,

t₀ is the time at which the peak of the wavepacket is located at a distance z₀ from the photon source, and

Δ² is the variance of the intensity spectrum.

The time to can be determined by z₀/v, where v is the velocity of the pulse traveling through free space or in an optical fiber.

The wavepacket ξ(t) is the amplitude of the pulse, and |ξ(t)|² is a photodetection probability density function of the pulse, where the photodetection probability density function |ξ(t)|² satisfies the normalization condition: ∫_(−∞)^(∞)𝕕tξ(t)² = 1 The probability of photodetection of a photon in the time interval (t₁,t₂) at a distance z₀ from the photon source is given by: Probability  of  (t₁ < t₂) = ∫_(t₁)^(t₂)𝕕tξ(t)²

FIG. 8 illustrates a probability distribution associated with a pulse output from a source 802 and transmitted in an optical fiber 804 to a detector 806. The horizontal line 808 represents the distance z₀ the photon travels from the source 802 to the detector 806, and the horizontal line 810 is a time axis. The curve 812 represents the photodetection probability density function |ξ(t)|². In FIG. 8, the photodetection probability density function |ξ(t)|² 812 is centered at the time to, which corresponds to the time a pulse takes to travel the distance z₀. An area under the curve 812 represents the probability of detecting the pulse within a particular time period. For example, hash-marked region 814 represents the probability of detecting the photon within the time period t₁<t₀<t₂. Time period 816 is called a “time bin” and corresponds to the time period within which the photon is detected at the detector 806.

The time dependent creation operators can be used to generate a photon wavepacket creation operator as follows: â_(ξ)^(†) = ∫_(−∞)^(∞)𝕕t  ξ(t)â^(†)(t) The creation operator can be used to construct continuous-mode number states that represent photons transmitted through free space or in an optical fiber as follows: ${\left. n_{\xi} \right\rangle = {\frac{\left( {\hat{a}}_{\xi}^{\dagger} \right)^{n}}{\sqrt{n!}}\left. 0 \right\rangle}},$ where |0

is the continuous-mode vacuum state. The continuous-mode number states satisfy the following same conditions: n̂n_(ξ)⟩ = nn_(ξ)⟩, ⟨n_(ξ)^(′)|n_(ξ)⟩ = δ_(n^(′)n), and ${\sum\limits_{n_{\xi} = 0}^{\infty}{\left. n_{\xi} \right\rangle\left\langle n_{\xi} \right.}} = 1$ As a result, the subscript ξ used to identify continuous-mode number states can be dropped. Note that the wavepacket constructed photon is not an eigenstate of any Hamiltonian.

An Overview of Coherent States

The most common kind of single-mode states are linear superpositions of the number states. There are a number of different possible linear superpositions of the number states, but the coherent state: $\left. \alpha \right\rangle = {{\exp\left( {{- \frac{1}{2}}{\alpha }^{2}} \right)}{\sum\limits_{n = 0}^{\infty}{\frac{\alpha^{n\quad}}{\sqrt{n!}}\left. n \right\rangle}}}$ is a linear superposition of the number states used in many applications of quantized electromagnetic radiation. The coherent states are eigenstates of the annihilation operator: â|α

=α|α

, where taking the complex conjugate gives:

α|â^(†)=

α|α* However, the coherent state |α

is not an eigenstate of the creation operator at because the summation over a cannot be rearranged to give the coherent state from at â^(†).

The coherent state expectation value for the number operator:

n

=

α|{circumflex over (n)}|α

=|α|² indicates that |α|² is the mean number of photons. The probability of detecting n photons in a measurement of the number of photons is a Poisson distribution: $P_{n} = {{\left\langle n \middle| \alpha \right\rangle }^{2} = {{\exp\left( {- {\alpha }^{2}} \right)}\frac{{\alpha }^{2}}{n!}}}$ The Poisson distribution approaches a Gaussian distribution for large values of |α|².

The coherent state is a quantum state whose properties most closely resemble a classical electromagnetic wave of stable amplitude and fixed phase. For example, the electric field operator corresponding to an electric field propagating in the z direction, with the mode subscripts k and s removed, is: $\begin{matrix} {{\hat{E}(\Omega)} = {\frac{1}{2}\left( {{\hat{a}{\mathbb{e}}^{- {\mathbb{i}\Omega}}} + {{\hat{a}}^{\dagger}{\mathbb{e}}^{\mathbb{i}\Omega}}} \right)}} \\ {= {{\hat{X}\cos\quad\Omega} + {\hat{Y}\sin\quad\Omega}}} \end{matrix}$ where the time t and displacement z are contained in the phase angle: ${{\Omega\left( {z,t} \right)} = {{\omega\quad t} - {kz} - \frac{\pi}{2}}},$ and the electric field is measured in units of √{square root over (

ω/2ε₀V)}.

The coherent state is a nearly classical state because it gives the correct sinusoidal form for the electric field expectation value or coherent signal: $\begin{matrix} {\Sigma = \left\langle {\alpha{{\hat{E}(\Omega)}}\alpha} \right\rangle} \\ {= {{\alpha }{\cos\left( {\Omega - \phi} \right)}}} \end{matrix}$ where α=|α|e^(iφ), and

φ is the mean phase angle of the coherent state excitation of the mode.

FIG. 9 is a phase-space diagram of a coherent state. Horizontal axis 902 is the X quadrature axis, and vertical axis 904 is the Y quadrature axis. Disk 906 represents an area of uncertainty associated with the phase and average number of photons of a coherent state |α

, the fluctuations being equal in all directions of phase space. Directional arrow 908 represents the amplitude |α|, or the mean number of photons, and (Ω−φ) is the angle above the X quadrature axis. The uncertainty disk is resolved into a photon-number uncertainty and phase uncertainty, as indicated by directional arrows 910 and 912, respectively. Directional arrow 910 represents uncertainty in the phase angle φ, which diminishes as |α| increases according to: ${\Delta\quad\phi} = \frac{1}{2{\alpha }}$ Because the phase uncertainty is inversely proportional to |α|, the coherent state becomes better defined as the average number of photons is increased. Directional arrow 912 represents the photon-number uncertainty: Δn=|α| The area of the uncertainty disk is: ${{\Delta\phi\Delta}\quad n} = \frac{1}{2}$ A measurement on the coherent state outputs the coherent signal Σ, which is represented by projecting the center of the uncertainty disk 906 onto the X-quadrature axis 902.

The angle Ω is a property of the measurement that can be set equal to zero by the experimentalist, which gives an X-quadrature-expectation value:

α|{circumflex over (X)}|α

=|α| cos φ=K A homodyne detection measurement outputs the X-quadrature-expectation value K as a function of the phase angle φ and the amplitude |α|. The quantity measured is the difference between the numbers of photons arriving at two different photodetectors during a period of time. The homodyne detector measures photon numbers, or photon counts, and the effect is to produce measurements proportional to an electric field quadrature, enabled through the measurement of detector currents.

A probability amplitude associated with homodyne detection of the coherent state is given by a Gaussian function: ${{f\left( {x,\beta} \right)} = {\frac{1}{\left( {2\pi} \right)^{1/4}}{\exp\left( {- \frac{\left( {x - {2\beta}} \right)^{2}}{4}} \right)}}},$ where β is equal to α cos φ, and

x is the value signal output from the homodyne detection.

Squaring the amplitude f(x,β) gives the probability distribution of a homodyne measurement. FIG. 10 is a plot of a probability distribution 1002 associated with the coherent state representation shown in FIG. 9. Horizontal axis 1004 corresponds to X quadrature values, and vertical axis 1006 is a probability axis. The probability distribution 1002 is centered at 2α cos φ. A coherent signal output by a homodyne measurement falls within the region 1008. The probability that the coherent signal falls within a particular interval is associated with an area under the probability distribution 1002. For example, hash-marked region 1010 represents the probability that the homodyne detection gives a value in the interval 1012.

An Overview of Quantum Entanglement

A quantum system comprising a first quantum subsystem and a second quantum subsystem has a Hilbert space H_(A)

H_(B), where H_(A) is a Hilbert space associated with the first quantum system, and H_(B) is a Hilbert space associated with the second quantum system. The kets |i

_(A) represent the orthonormal eigenstates of the Hilbert space H_(A), and the kets |j

_(B) represents the orthonormal eigenstates of the Hilbert space H_(B), where i and j are positive integers. Any linear superposition of states in the Hilbert space H_(A)

H_(B) is given by: ${\left. \Psi \right\rangle_{AB} = {\sum\limits_{i,j}{c_{ij}\left. i \right\rangle_{A}\left. j \right\rangle_{B}}}},$ where the amplitudes c_(ij) are complex numbers satisfying the condition: ${\sum\limits_{ij}{c_{ij}}^{2}} = 1$

Special kinds of linear superpositions of states |Ψ

_(AB) are called “direct product states” and are represented by the product: ${\left. \Psi \right\rangle_{AB} = {{\left. \psi \right\rangle_{A}\left. \psi \right\rangle_{B}} = {\left( {\sum\limits_{i}{c_{i}^{(A)}\left. i \right\rangle_{A}}} \right)\left( {\sum\limits_{j}{c_{j}^{(B)}\left. j \right\rangle_{B}}} \right)}}},$ where

|ψ

_(A) is a normalized linear superposition of states in the Hilbert space H_(A);

and |ψ

_(B) is a normalized linear superposition of states in the Hilbert space H_(B). However, linear superpositions in the Hilbert space H_(A)

H_(B) that cannot be written as a product state are entangled states. In general, for a Hilbert space comprising two or more quantum subsystems, an entangled state is a linear superposition of states that cannot be written as a direct product state. The set of GHZ states are examples of entangled qubit states. For example, consider an entangled 4-qubit, GHZ state: $\left. \psi \right\rangle_{4004} = {\frac{1}{\sqrt{2}}\left( {{\left. 0 \right\rangle_{1}\left. 0 \right\rangle_{2}\left. 0 \right\rangle_{3}\left. 0 \right\rangle_{4}} + {\left. 1 \right\rangle_{1}\left. 1 \right\rangle_{2}\left. 1 \right\rangle_{3}\left. 1 \right\rangle_{4}}} \right)}$ Suppose a measurement performed on the first qubit system in the basis {|0

₁, |1

₁} projects the state of the first qubit system onto the basis qubit state |0

₁. According to the state |ψ

₄₀₀₄, the measurement also projects the state of the second, third, and fourth qubit systems onto the basis qubit states |0

₂, |0

₃, and |0

₄, respectively. In addition, the state |ψ

₄₀₀₄ cannot be written as a product of the qubits α₁|0

₁+β₁|1

₁, α₂0

₂+β₂|1

₂, α₃|0

₃+β₃|1

₃, and α₄|0₄

₄+β₄|1

₄, for any choice of the parameters α₁, β₁, α₂, β₂, α₃, β₃, α₄, and β₄.

Embodiments of the Present Invention

Various embodiments of the present invention are directed to methods for determining a phase shift φ acquired by entangled N-qubit systems. The entangled N-qubit systems are represented by a GHZ NOON state: $\left. \psi \right\rangle_{NOON} = {\frac{1}{\sqrt{2}}\left( {\left. {N,0} \right\rangle + \left. {0,N} \right\rangle} \right)}$ where

|N, 0

represents the product state |0

₁ . . . |0

_(N); and

|0,N

represents the product state |1

₁ . . . |1

_(N). The entangled N-qubit systems accumulate a relative phase e^(iNφ) when each qubit system in the state |1

acquires the phase shift φ, which is represented by the state: $\left. \psi_{N} \right\rangle = {\frac{1}{\sqrt{2}}\left( {\left. {N,0} \right\rangle + {{\mathbb{e}}^{{\mathbb{i}}\quad N\quad\varphi}\left. {0,N} \right\rangle}} \right)}$ The methods of the present invention separately interact each of the N-qubit systems with a probe electromagnetic field, but the N-qubit systems are kept separate to prevent qubit-system interactions. In different embodiments of the present invention, the qubit interactions with the probe electromagnetic field causes the probe electromagnetic field to acquire either a phase shift or a change in the average number of photons, which can then used to determine the phase shift φ. For the sake of simplicity, methods of the present invention are described below with reference to an entangled 4-qubit system represented by the state: $\left. \psi_{4} \right\rangle = {\frac{1}{\sqrt{2}}\left( {\left. 0000 \right\rangle + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle}} \right)}$ where the relative phase e^(i4φ) is accumulated when each qubit system in the state |1

of the tensor product |1111

acquires the unknown and sought after phase shift φ.

In a first embodiment of the present invention, the state |ψ₄

represents four entangled matter-based qubit systems that are each represented separately by the qubits |ψ

₁, |ψ

₂, |ψ

₃, and |ψ

₄. The qubits can represent, for example, two discrete energy levels of an atom or the spin-up and spin-down states of an electron. FIG. 11 illustrates a first implementation 1100 for determining the phase shift φ acquired by the state |ψ₄

that represents an embodiment of the present invention. The implementation 1100 comprises an interaction region 1102, homodyne detection 1104, four separate Hadamard gates 1106-1109, and four qubit detectors 1112-1115. The interaction region 1102 comprises four separate and identical interaction mediums 1120-1123 described below. Transmission channels 1126-1129 separately transmit four qubit systems to the interaction mediums 1120-1123, respectively, and to the Hadamard gates 1106-1109, respectively. The separate transmission channels 1126-1129 prevent the qubit systems from interacting with one another. A transmission channel 1130 transmits an electromagnetic field represented by a coherent state |α

to a beamsplitter 1132, which divides the electromagnetic field into a probe electromagnetic field and a reference electromagnetic field. The probe electromagnetic field is represented by a coherent state |α

_(a) of |α|² and is called a “probe field.” The probe field is transmitted to the interaction region 1102 and separately interacts with each of the qubit systems at each of the interaction mediums 1120-1123, respectively. The reference electromagnetic field is represented by a coherent state |α

_(a) ^(ref) and is called a “reference field.” The reference field is directed by mirrors 1134 and 1136 to the homodyne detection 1104. The entangled 4-qubit systems and the probe field input to the interaction region 1102 are represented by an input state: $\begin{matrix} {\left. \Psi \right\rangle_{i\quad n} = {\left. \psi_{4} \right\rangle\left. \alpha \right\rangle_{a}}} \\ {= {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. \alpha \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. \alpha \right\rangle_{a}}} \right)}} \end{matrix}$

The probe field and the reference field both have identical average photon numbers and identical phases. The probe field acquires a phase shift, represented by g(θ), as a result of interacting separately with each of the qubit systems at the interaction mediums 1120-1123. A resultant phase-shifted probe field, represented by the state |αe^(ig(θ))

_(a), is output from the interaction region 1102. The phase shift g(θ) can be determined by comparing the state |αe^(ig(θ))

_(a) of the phase-shifted probe field to the state |α

_(a) ^(ref) of the reference field at the homodyne detection 1104. As a result, the probe field input to the interaction region 1102 is initially thought of as having a phase angle equal to “0.”

FIG. 12 is a phase space diagram of the probe field input to the interaction region 1102, in FIG. 11, that represents an embodiment of the present invention. In FIG. 12, horizontal axis 1202 represents the continuous spectrum of X-quadrature expectation values, vertical axis 1204 represents the continuous spectrum of Y-quadrature expectation values, directional arrow 1206 represents the probe field, and disk 1208 represents the uncertainty in the phase angle and the uncertainty in the average number of photons in the state |α

_(a) of the probe field. Note that the amplitude |α|, or length, of the directional arrow 1206 represents an average number of photons in the probe field state |α

_(a).

In FIG. 11, the interaction mediums 1120-1123 couple the probe field and the corresponding qubit system. The interaction mediums 1120-1123 each produce an electric, or magnetic, dipole interaction between the probe field and the corresponding qubit system that is represented by an interaction Hamiltonian: Ĥ _(j) ⁽¹⁾=

χ(â ^(†) â){circumflex over (σ)}_(z,j) where

(â^(\)â) is a number operator that operates on the number states of the probe field;

χ is a constant that represents the coupling strength between the probe field and the qubit system;

superscript “1” refers to the first embodiment;

subscript j is a qubit index; and

{circumflex over (σ)}_(z,j)=|0

_(j j)

0|−|1 is an inversion operator.

The inversion operator {circumflex over (σ)}_(z,j) operates on the qubit basis states as follows: {circumflex over (σ)}_(z,j)|0

_(i)=δ_(ij)|0

_(i), and {circumflex over (σ)}_(z,j)|1

_(i)=−δ_(ij)|1

_(i)

The dipole interaction at each of the interaction mediums 1120-1123 occurs for a period of time t. Each dipole interaction is mathematically represented by a dipole time-evolution operator: $\begin{matrix} {{\hat{U}}_{j}^{(1)} = {\exp\left( {{- {\mathbb{i}}}\quad{\hat{H}}_{j}^{1}{t/\hslash}} \right)}} \\ {= {\exp\left( {{- {{\mathbb{i}\theta}\left( {{\hat{a}}^{\dagger}\hat{a}} \right)}}{\hat{\sigma}}_{z,j}} \right)}} \end{matrix}$ where

θis equal to the product χt.

The qubit systems pass unchanged through each of the interaction mediums 1120-1123. However, the dipole interaction between a single qubit system and the probe field, at an interaction medium, causes the probe field to acquire either a positive or a negative phase shift θ depending on whether the corresponding qubit system transmitted to the interaction medium is in the state |0

or in the state |1

. For example, at the interaction medium 1120, a dipole interaction between a qubit system in the state |0

₁ and the probe field in the state |α

_(a) produces a phase-shifted probe field is the state |αe^(iθ)

_(a), which is mathematically represented by applying the dipole time-evolution operator Û₁ ⁻¹ as follows: Û₁ ⁽¹⁾|0

₁|α

_(a)=|0

₁|αe^(iθ)

_(a)

The dipole interaction between the qubit system in the state |1

₁ and the probe field in the state |α

_(a) at the interaction medium 1120 produces a phase-shifted probe field in the state |αe^(−iθ)

_(a), which is also mathematically represented by: Û₁ ⁽¹⁾|1

₁|α

_(a)=|1

₁|αe^(−iθ)

_(a) Note that the dipole time-evolution operator Û₁ ⁻¹ does not change the state of the qubit basis states |0

and |1

.

The interaction region 1102 receives a quantum system represented by the input state |Ψ

_(in) and outputs an entangled probe-qubit system represented by an entangled output state |Ψ

_(out). The output state |Ψ

_(out) is determined by successively applying the dipole time-evolution operator corresponding to each dipole interaction at the interaction mediums 1120-1123 as follows. First, the dipole interaction at the interaction medium 1120 produces a first entangled probe-qubit system whose state is determined as follows: $\begin{matrix} {{{\hat{U}}_{1}^{(1)}\left. \Psi \right\rangle_{i\quad n}} = {\frac{1}{\sqrt{2}}\left( {{{\hat{U}}_{1}^{(1)}\left. 0000 \right\rangle\left. \alpha \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}{\hat{U}}_{1}^{(1)}\left. 1111 \right\rangle\left. \alpha \right\rangle_{a}}} \right)}} \\ {= {\frac{1}{\sqrt{2}}\begin{pmatrix} {{\left( {{\hat{U}}_{1}^{(1)}\left. 0 \right\rangle_{1}\left. \alpha \right\rangle_{a}} \right)\left. 0 \right\rangle_{2}\left. 0 \right\rangle_{3}\left. 0 \right\rangle_{4}} +} \\ {{{\mathbb{e}}^{\mathbb{i}4\varphi}\left( {{\hat{U}}_{1}^{(1)}\left. 1 \right\rangle_{1}\left. \alpha \right\rangle_{a}} \right)}\left. 1 \right\rangle_{2}\left. 1 \right\rangle_{3}\left. 1 \right\rangle_{4}} \end{pmatrix}}} \\ {= {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha\mathbb{e}}^{\mathbb{i}\theta} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha\mathbb{e}}^{- {\mathbb{i}\theta}} \right\rangle_{a}}} \right)}} \end{matrix}$ Next, in a similar manner, the first entangled probe-qubit system is input to the interaction medium 1121 produces a second entangled probe-qubit system whose state is represented by: ${{\hat{U}}_{2}^{(1)}{\hat{U}}_{1}^{(1)}\left. \Psi \right\rangle_{i\quad n}} = {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha\mathbb{e}}^{\mathbb{i}2\theta} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha\mathbb{e}}^{- {\mathbb{i}2\theta}} \right\rangle_{a}}} \right)}$ Next, the second entangled probe-qubit system is input to the interaction medium 1122, which produces a third entangled probe-qubit quantum system whose state is represented by: ${{\hat{U}}_{3}^{(1)}{\hat{U}}_{2}^{(1)}{\hat{U}}_{1}^{(1)}\left. \Psi \right\rangle_{i\quad n}} = {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha\mathbb{e}}^{\mathbb{i}3\theta} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha\mathbb{e}}^{- {\mathbb{i}3\theta}} \right\rangle_{a}}} \right)}$ Finally, the third entangled probe-qubit quantum system is input to the interaction medium 1123, which produces a fourth entangled probe-qubit system whose state is represented by: $\begin{matrix} {\left. \Psi \right\rangle_{out} = {{\hat{U}}_{4}^{(1)}{\hat{U}}_{3}^{(1)}{\hat{U}}_{2}^{(1)}{\hat{U}}_{1}^{(1)}\left. \Psi \right\rangle_{i\quad n}}} \\ {= {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha\mathbb{e}}^{\mathbb{i}4\theta} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha\mathbb{e}}^{- {\mathbb{i}4\theta}} \right\rangle_{a}}} \right)}} \end{matrix}$

FIG. 13 shows four successive and separate phase-space diagrams 1302-1305 associated with a positive phase shift 9 acquired by the probe field at each of the interaction mediums 1120-1123, in FIG. 11, that represent an embodiment of the present invention. Each diagram, in FIG. 13, includes a dashed-line disk, such as disk 1306 in diagram 1302, that represents the uncertainty disk 1208 associated with the probe field initially in the state |α

_(a), in FIG. 12. Diagrams 1302-1305 shows successive phase changes in the state of the probe field as a result of each dipole interaction between the probe field and the qubit systems in the tensor product state |0000

at the interaction mediums 1120-1123, respectively, in FIG. 11. For example, the probe field input to the interaction medium 1120 is initially in the state |α_(a) and is output from the interaction medium 1120 in the phase-shifted state |αe^(iθ)

_(a) as shown in diagram 1302.

FIG. 14 shows four successive and separate phase-space diagrams 1402-1405 associated with a negative phase shift θ acquired by the probe field at each interaction medium, in FIG. 11, that represent embodiments of the present invention. Each diagram, in FIG. 14, includes a dashed-line disk, such as disk 1406 in diagram 1402, that represents the uncertainty disk 1208 associated with the probe field initially in the state |αa

_(a), in FIG. 12. Diagrams 1402-1405 shows successive phase changes in the state of the probe field as a result of each dipole interaction between the probe field and the qubit systems in the tensor product state |1111

at the interaction mediums 1120-1123, respectively, in FIG. 11. For example, the probe field input to the interaction medium 1120 is initially in the state |α

_(a) and is output from the interaction medium 1120 in the phase-shifted state |αe^(−iθ)

_(a) as shown in diagram 1402.

The interaction region 1102, in FIG. 11, outputs an entangled probe-qubit system represented by the output state: $\left. \Psi \right\rangle_{out} = {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha\mathbb{e}}^{\mathbb{i}4\theta} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha\mathbb{e}}^{- {\mathbb{i}4\theta}} \right\rangle_{a}}} \right)}$

The qubit systems are input to the Hadamard gates 1106-1109. Each of the Hadamard gates 1106-1109 performs a qubit system operation that is mathematically represented by the Hadamard operator: ${\hat{H}}_{Had} = {\frac{1}{\sqrt{2}}\left( {{\left. 0 \right\rangle\left\langle 0 \right.} + {\left. 1 \right\rangle\left\langle 0 \right.} + {\left. 0 \right\rangle\left\langle 1 \right.} - {\left. 1 \right\rangle\left\langle 1 \right.}} \right)}$ The Hadamard gate operator operates on the qubit basis states |0

and |1

as follows: ${{{\hat{H}}_{Had}\left. 0 \right\rangle} = {\frac{1}{\sqrt{2}}\left( {\left. 0 \right\rangle + \left\langle 1 \right.} \right)}},{and}$ ${{\hat{H}}_{Had}\left. 1 \right\rangle} = {\frac{1}{\sqrt{2}}\left( {\left. 0 \right\rangle - \left\langle 1 \right.} \right)}$

The Hadamard gates 1106-1109 are used, with subsequent qubit measurements, to disentangle the probe field in the phase-shifted probe field states |αe^(i4θ)

and |αe^(i4θ)

_(a) from the qubit systems and transfer the phase shift φ to a linear superposition of the phase-shifted probe field states |αe^(−i4θ)

_(a) and |e^(i4θ)

_(a). A Hadamard output state, |Ψ

_(Had), represents the quantum system output from the Hadamard gates 1106-1109 and is obtained as follows: $\begin{matrix} {\left. \Psi \right\rangle_{Had} = {{\hat{H}}_{Had}^{\otimes 4}\left. \Psi \right\rangle_{out}}} \\ {= {\frac{1}{\sqrt{2}}\left( {{{\hat{H}}_{Had}^{\otimes 4}\left. 0000 \right\rangle\left. {\alpha\mathbb{e}}^{\mathbb{i}4\varphi} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}{\hat{H}}_{Had}^{\otimes 4}\left. 1111 \right\rangle\left. {\alpha\mathbb{e}}^{\mathbb{i}4\theta} \right\rangle_{a}}} \right)}} \end{matrix}$

Expanding the Hadamard gate operations Ĥ_(Had) ⁴|0000

and Ĥ_(Had) ⁴|0000

and rearranging the product states gives the Hadamard output state: $\left. \Psi \right\rangle_{Had} = {{\left( {{\frac{1}{4}\left. \Gamma \right\rangle_{1}} + {\frac{\sqrt{6}}{4}\left. \Gamma \right\rangle_{3}} + {\frac{1}{4}\left. \Gamma \right\rangle_{5}}} \right)\left. \xi \right\rangle_{field}^{+}} + {\left( {{\frac{1}{2}\left. \Gamma \right\rangle_{2}} + {\frac{1}{2}\left. \Gamma \right\rangle_{4}}} \right)\left. \xi \right\rangle_{field}^{-}}}$ where ${\left. \Gamma \right\rangle_{1} = \left. 0000 \right\rangle},{\left. \Gamma \right\rangle_{2} = {\frac{1}{2}\left( {\left. 1000 \right\rangle + \left. 0100 \right\rangle + \left. 0010 \right\rangle + \left. 0001 \right\rangle} \right)}},{\left. \Gamma \right\rangle_{3} = {\frac{1}{\sqrt{6}}\left( {\left. 1100 \right\rangle + \left. 1010 \right\rangle + \left. 1001 \right\rangle + \left. 0110 \right\rangle + \left. 0101 \right\rangle + \left. 0011 \right\rangle} \right)}},{\left. \Gamma \right\rangle_{4} = {\frac{1}{2}\left( {\left. 1110 \right\rangle + \left. 1101 \right\rangle + \left. 1011 \right\rangle + \left. 0111 \right\rangle} \right)}},{\left. \Gamma \right\rangle_{5} = \left. 1111 \right\rangle},{\left. \xi \right\rangle_{field}^{1 +} = {\frac{1}{\sqrt{2}}\left( {\left. {\alpha\mathbb{e}}^{\mathbb{i}4\theta} \right\rangle + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {\alpha\mathbb{e}}^{- {\mathbb{i}4\theta}} \right\rangle}} \right)}},{and}$ $\left. \xi \right\rangle_{field}^{1 -} = {\frac{1}{\sqrt{2}}\left( {\left. {\alpha\mathbb{e}}^{\mathbb{i}4\theta} \right\rangle - {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {\alpha\mathbb{e}}^{- {\mathbb{i}4\theta}} \right\rangle}} \right)}$ The states |ξ

_(field) ¹⁺ and |ξ

_(field) ¹⁻ are called “field states,” and each field state represents the probe field in a different coherent linear superposition of the states |αe^(−i4θ)

_(a) and |αe^(i4θ)

_(a). Note that as a result of the Hadamard gate operations on the qubit systems in the output state |Ψ

_(out), the field states |ξ

_(field) ¹⁺ and |ξ

_(field) ¹⁻ have acquired the phase shift φ.

Whether the probe field is in the field state |ξ

_(field) ¹⁺ or in the field state |ξ

_(field) ¹⁻ when it reaches the homodyne detection 1104 is entirely random. However, the field state of the probe field that reaches the homodyne detection 1104 is determined by performing qubit measurements at the qubit detectors 1112-1115. The qubit measurements project the state of the probe field onto either the field state |ξ

_(field) ¹⁺ or the field state |ξ

_(field) ¹⁻. For example, suppose the quantum system represented by the Hadamard output state |Ψ

_(Had) reaches the qubit detectors 1112-1115, and the qubit measurements reveal that the qubits are in the tensor product state |0100

of the state |Γ

₂ above. According to the Hadamard output state |Ψ

_(Had), the qubit measurements project the Hadamard output state |Ψ

_(Had) onto the state |1010

|ξ

_(field) ¹⁺, which means that the state of the probe field just after the qubit measurements is the field state |ξ

_(field) ¹⁺. The homodyne detection 1104 measures the electric field intensity, or X-quadrature-expectation value, of the probe field as functions of the phase angle associated with the phase-shifted probe fields |αe^(−i4θ)

_(a) and |αe^(i4θ)

_(a).

In order to determine the phase shift p acquired by the probe field in the field states |ξ

_(field) ¹⁺ and |ξ

_(field) ¹⁻, a large number of input states |Ψ

_(in) are input to the implementation 1100. After each of the qubit measurements, the homodyne detection 1104 X-quadrature-expectation values are stored in two of different data sets according to whether the probe field is in the field state |ξ

_(field) ¹⁺ or in the field state |ξ

_(field) ¹⁻. For example, when a qubit measurement projects the state of the probe field onto the field state |ξ

_(field) ¹⁺, the corresponding X-quadrature-expectation value is stored in a first data set, and when a subsequent qubit measurement projects the state of the probe field onto the field state |ξ

_(field) ¹⁻, the corresponding X-quadrature-expectation value is stored in a second data set. The homodyne detection 1104 X-quadrature-expectation values can be stored electronically, magnetically, or optically and can be processed using logical operations encoded in computer programs.

The first data set and the second data set are used to produce a first frequency distribution and a second frequency distribution, respectively. FIGS. 15A-15B are diagrams of example frequency distributions associated with a large number of hypothetical homodyne-detection results obtained for the probe field in the field states |ξ

_(field) ¹⁺ and |ξ

_(field) ¹⁻, respectively, that represent embodiments of the present invention. The example frequency distributions are obtained from a large number of hypothetical executions of the implementation 1100 described above with reference to FIG. 11. In FIGS. 15A-15B, horizontal axes, such as horizontal axis 1502 in FIG. 15A, represent continuous spectrums of X-quadrature-expectation values, and vertical axes, such as vertical axis 1504 in FIG. 15A, are normalized number of incidents axes. Each point represents a fraction of X-quadrature-expectation values that fall within a particular X-quadrature-expectation-value interval. The value associated with each point can be determined by averaging the number of the X-quadrature-expectation values that fall within each of the X-quadrature-expectation-value intervals. For example, point 1506 represents the fraction of X-quadrature-expectation values that fall within X-quadrature-expectation-value interval 1508. After a large number of hypothetical executions of the implementation 1100, the fraction associated with each point can also be thought of as the probability of a subsequent X-quadrature-expectation value falling within a particular X-quadrature-expectation-value interval. For example, the fraction associated with the point 1506 is the probability that a subsequent X-quadrature-expectation value associated with a probe field in the field state |ξ

_(field) ¹⁺ falls within the X-quadrature-expectation-value interval 1508.

Either of the frequency distributions shown in FIGS. 15A-15B can be used to determine the phase shift φ. Consider using the first frequency distribution in FIG. 15A. Applying the homodyne-detection operator |x

x| to the field state |ξ

_(field) ¹⁺ gives: $\begin{matrix} {{\left. x \right\rangle\left\langle x \middle| \xi \right\rangle_{field}^{+}} = {{F^{+}(x)}\left. x \right\rangle}} \\ {{= {\frac{1}{\sqrt{2}}{f(x)}\left( {{\mathbb{e}}^{{\mathbb{i}\phi}{(x)}} + {{\mathbb{e}}^{\mathbb{i}4\phi}{\mathbb{e}}^{- {{\mathbb{i}\phi}{(x)}}}}} \right)\left. x \right\rangle}},} \end{matrix}$ where ${{f(x)} = {\frac{1}{\left( {2\pi} \right)^{1/4}}{\exp\left( {- \frac{\left( {x - {2\alpha\quad\cos\quad 4\quad\theta}} \right)^{2}}{4}} \right)}}},{and}$ φ(x)=α sin 4θ(x−2α cos 4θ)mod2π.

The coefficient F⁺(x) is called a “probability amplitude,” which can be used to determine the probability distribution associated with the field state |ξ

_(field) ¹⁺ as follows: $\begin{matrix} {{P^{+}(x)} = {{F^{+}(x)} \cdot {F^{+ *}(x)}}} \\ {= {{f^{2}(x)}\left( {1 + {\cos\left( {{2{\phi(x)}} - {4\varphi}} \right)}} \right)}} \end{matrix}$ A least squares method, or other data fitting methods, can be used to fit the probability distribution P⁺(x) to the corresponding frequency distribution, in FIG. 15A, as function of the phase shift q). The probability distribution P⁺(x) is least squares fit to the points in the first frequency distribution by minimizing the sum of the squares of the distances between the frequency-distribution points and the probability distribution P⁺(x) as a function of the phase shift φ. The optimal phase shift φ determined by the least squares method gives a best fit probability distribution P⁺(x) to the frequency-distribution points. The least squares method, or other data fitting methods, can be used to determine the optimal phase shift φ using the probability distribution P⁻(x) and the second frequency distribution, in FIG. 15B.

FIG. 16 is a diagram of the probability distribution P⁺(x) fit to the frequency distribution shown in FIG. 15A that represents an embodiment of the present invention. In FIG. 16, curve 1602 represents the probability distribution P⁺(x) fit to the points of the frequency distribution in FIG. 15A for a hypothetical optimal phase shift φ.

In a second embodiment of the present invention, the state |ψ₄

represents four entangled optical qubit systems, each represented separately by the qubits |ψ

₁, |ψ

₂, |ψ

₃, and |ψ

₄. An optical qubit system can, for example, represent the two photon states of an electromagnetic field or vertically and horizontally polarized photons. FIG. 17 illustrates a second implementation 1700 for determining the phase shift φ acquired by the state |ψ₄

that represents an embodiment of the present invention. The implementation 1700 comprises an interaction region 1702, rotate quadrature operation 1704, homodyne detection 1706, four separate Hadamard gates 1708-1711, and four qubit detectors 1712-1715. The interaction region 1702 comprises four separate and identical cross-Kerr interaction mediums 1720-1723 described below. Transmission channels 1726-1729 separately transmit the four qubit systems to the interaction mediums 1720-1723 and to the Hadamard gates 1708-1711, respectively. The separate transmission channels 1726-1729 prevent the qubit systems from interacting with one another. A transmission channel 1730 receives an electromagnetic field represented by a coherent state |α

. A beamsplitter 1732 divides the electromagnetic field into a probe field and a reference field as described above with reference to FIG. 11. The probe field is represented by coherent state |α

_(a), and the reference field is represented by the coherent state |α

_(a) ^(ref). The reference field is directed by mirrors 1734 and 1736 to the homodyne detection 1706. The four qubit systems and the probe field are input to the interaction region 1702 and are represented by the input state |Ψ

_(in) described above with reference to FIG. 11. The probe field interacts separately with each of the qubit systems at each of the interaction mediums 1720-1723, respectively. As described above with reference to FIG. 12, the probe field input to the interaction region 1702 is initially thought of as having a phase angle equal to “0.”

The interaction mediums 1720-1723 separately couple the probe field to each qubit system. The interaction at each of the cross-Kerr interaction mediums 1720-1723 is represented by a cross-Kerr interaction Hamiltonian: Ĥ _(j) ⁽²⁾=

χ(â^(†)â)({circumflex over (b)}^(†){circumflex over (b)})_(j) where

superscript “2” refers to the second embodiment; and

({circumflex over (b)}^(†){circumflex over (b)}) is a number operator that operates on the qubits.

The number operator ({circumflex over (b)}^(†){circumflex over (b)})_(j) operates on the photon basis states as follows: ({circumflex over (b)}^(†){circumflex over (b)})_(j)|0

_(j)=0, and ({circumflex over (b)}^(†){circumflex over (b)})_(j)|1

_(j)=1·|1

_(j)

The cross-Kerr interaction at each of the interaction mediums 1720-1723 occurs for a period of time t and is mathematically represented by a cross-Kerr time-evolution operator: $\begin{matrix} {{\hat{U}}_{j}^{(2)} = {\exp\left( {{- {\mathbb{i}}}\quad{\hat{H}}_{j}^{2}{t/\hslash}} \right)}} \\ {= {\exp\left( {{- {{\mathbb{i}\theta}\left( {{\hat{a}}^{\dagger}\hat{a}} \right)}}\left( {{\hat{b}}^{\dagger}\hat{b}} \right)_{j}} \right)}} \end{matrix}$ where

θ is equal to the product χt.

The qubit systems pass unchanged through each of the interaction mediums 1720-1723. However, the cross-Kerr interaction between a qubit system in the state |0

and the probe field leaves both the qubit |0

and the probe field in the coherent state |α

_(a) unchanged, but the cross-Kerr interaction between a qubit system in the state |1

and the probe field causes the probe field to acquire a negative phase shift θ. For example, at the interaction medium 1120, a cross-Kerr interaction between a qubit system in the state |1

₁ and the probe field in the coherent state |α

_(a) produces a phase-shifted probe field state |αe^(−iθ)

_(a) which is mathematically determined by: Û ₁ ⁽²⁾|1

₁|α

_(a)=|1

₁ |αe ^(−iθ)

_(a) Note that the cross-Kerr time-evolution operator Û_(j) ² does not change the state of the qubit basis states.

The interaction region 1702 outputs an entangled probe-qubit system represented by the output state: $\left. \Psi \right\rangle_{out} = {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. \alpha \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha\mathbb{e}}^{- {\mathbb{i}4\theta}} \right\rangle_{a}}} \right)}$ The output state |Ψ

_(out) can be determined mathematically by successively and separately applying the cross-Kerr time-evolution operators associated with each of the cross-Kerr interaction mediums 1720-1723, as described above with reference to the interaction region 1102 in FIG. 11.

The probe fields |α

_(a) and |αe^(−i4θ)

_(a) in the output state |Ψ

_(out) are transmitted to the rotate quadrature operation 1704, which produces a probe field phase shift of 2θ. The rotate quadrature operation 1704 is mathematically represented by a rotation operator, exp (i2θâ^(†)â), which when applied to the output state |Ψ

_(out) gives the state: $\begin{matrix} {\left. \Psi \right\rangle_{out}^{rot} = {{\exp\left( {{\mathbb{i}}\quad 2\theta\quad{\hat{a}}^{\dagger}\hat{a}} \right)}\left. \Psi \right\rangle_{out}}} \\ {= {\frac{1}{\sqrt{2}}\begin{pmatrix} {{\left. 0000 \right\rangle{\exp\left( {{\mathbb{i}}\quad 2\quad\theta\quad{\hat{a}}^{\dagger}\hat{a}} \right)}\left. \alpha \right\rangle_{a}} +} \\ {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle{\exp\left( {{\mathbb{i}}\quad 2\theta\quad{\hat{a}}^{\dagger}\hat{a}} \right)}\left. {\alpha\mathbb{e}}^{- {\mathbb{i}4\theta}} \right\rangle_{a}} \end{pmatrix}}} \\ {= {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha\mathbb{e}}^{\mathbb{i}2\theta} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha\mathbb{e}}^{- {\mathbb{i}2\theta}} \right\rangle_{a}}} \right)}} \end{matrix}$

FIGS. 18A-18B are phase-space diagrams of 2θ phase-shifted probe fields output from the rotate quadrature operation 1704, in FIG. 17, that represents an embodiment of the present invention. In FIG. 18A, the probe field, represented by a dashed-line directional arrow 1802 and uncertainty disk 1804, is transformed into phase-shifted probe field state |αe^(i2θ)

_(a), which is represented by directional arrow 1806 and uncertainty disk 1808. In FIG. 18B, the phase-shifted probe field state |αe^(−i4θ)

_(a), represented by a dashed-line directional arrow 1810 and uncertainty disk 1812, is transformed into a phase-shifted probe field state |αe^(i2θ)

_(a), which is represented by directional arrow 1814 and uncertainty disk 1816.

The qubits represented in the output state |Ψ

_(out) ^(rot) are input to the Hadamard gates 1708-1711. The Hadamard gates 1708-1711, with subsequent qubit measurements, disentangle the phase-shifted probe fields from the entangled qubit systems represented by the output state |Ψ

_(out) ^(rot). A Hadamard output state output from the Hadamard gates 1708-1711 is mathematically determined as described above with reference to FIG. 11 and is given by: $\begin{matrix} {\left. \Psi \right\rangle_{Had} = {{\hat{H}}_{Had}^{\otimes 4}\left. \Psi \right\rangle_{out}^{rot}}} \\ {= {{\left( {{\frac{1}{4}\left. \Gamma \right\rangle_{1}} + {\frac{\sqrt{6}}{4}\left. \Gamma \right\rangle_{3}} + {\frac{1}{4}\left. \Gamma \right\rangle_{5}}} \right)\left. \xi \right\rangle_{field}^{2 +}} + {\begin{pmatrix} {{\frac{1}{2}\left. \Gamma \right\rangle_{2}} +} \\ {\frac{1}{2}\left. \Gamma \right\rangle_{4}} \end{pmatrix}\left. \xi \right\rangle_{field}^{2 -}}}} \end{matrix}$ where ${\left. \xi \right\rangle_{field}^{2 +} = {\frac{1}{\sqrt{2}}\left( {\left. {\alpha\mathbb{e}}^{\mathbb{i}2\theta} \right\rangle + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {\alpha\mathbb{e}}^{- {\mathbb{i}2\theta}} \right\rangle}} \right)}},{and}$ $\left. \xi \right\rangle_{field}^{2 -} = {\frac{1}{\sqrt{2}}{\left( {\left. {\alpha\mathbb{e}}^{\mathbb{i}2\theta} \right\rangle - {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {\alpha\mathbb{e}}^{- {\mathbb{i}2\theta}} \right\rangle}} \right).}}$

The field states |ξ

_(field) ²⁺ and |ξ

_(field) ²⁻ represent two coherent linear superpositions of the probe field states, which have both acquired the phase shift (P. The field state that reaches the homodyne detection 1104 is determined by the qubit measurements performed at the qubit detectors 1712-1715. The qubit measurements project the state of the probe field onto either the field state |ξ

_(field) ²⁺ or onto the field state |ξ

_(field) ²⁻, as described above with reference to FIG. 11.

In order to determine the phase shift φ acquired by the probe field in the states |ξ

_(field) ²⁺ and |ξ

_(field) ²⁻, a large number of identically prepared input states |Ψ

_(in) are successively input to the implementation 1700. After each of the qubit measurements, X-quadrature-expectation values collected from the homodyne detection 1706 are stored in a first data set when the probe field is found in the field state |ξ

_(field) ²⁺, and X-quadrature-expectation values are stored in a second data set when the probe field is found in the field state |ξ

_(field) ²⁻. The first data set is used to produce a first frequency distribution associated with the field state |ξ

_(field) ²⁺, and the second data set is used to produce a second frequency distribution associated with the field state |ξ

_(field) ²⁻, as described above with reference to FIGS. 15A-15B.

As described above with reference to FIGS. 15A-16, either the first or the second frequency distribution can be used to determine the phase shift φ. For example, the probability amplitudes associated with the probe field in the field state |ξ

_(field) ²⁺ is: ${F^{+}(x)} = {\frac{1}{\sqrt{2}}{f(x)}\left( {{\mathbb{e}}^{{\mathbb{i}\phi}{(x)}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}{\mathbb{e}}^{- {{\mathbb{i}\phi}{(x)}}}}} \right)}$ where ${{f(x)} = {\frac{1}{\left( {2\pi} \right)^{1/4}}{\exp\left( {- \frac{\left( {x - {2\alpha\quad\cos\quad 2\theta}} \right)^{2}}{4}} \right)}}},{and}$ ϕ(x) = α  sin   2θ(x − 2α  cos   2θ)  mod  2π. The probability distribution associated with the probe field in the field state |ξ

_(field) ²⁺ is: P ⁺(x)=f ²(x)(1+cos(2φ(x)−4φ)).

An optimal phase shift φ value can be determined by applying a least squares method, or other curve fitting method, to fit the probability distribution P⁺(x) to the data points of the first frequency distribution, as describe above with reference to FIG. 16. Note also the probability distribution P⁻(x) associated with the second frequency distribution can also be used to determine the optical phase shift (value.

FIG. 19 illustrates a third implementation 1900 for determining the phase shift φ acquired by qubit systems represented by the state |ψ₄

that represents an embodiment of the present invention. The implementation 1900 comprises an interaction region 1902, displacement operation 1904, photon number counter 1906, four separate Hadamard gates 1908-1911, and four qubit detectors 1912-1915. The interaction region 1902 comprises four separate and identical interaction mediums 1920-1923 described below. Transmission channels 1926-1929 separately transmit four qubit systems to the interaction mediums 1920-1923, respectively. The separate transmission channels 1926-1929 prevent the qubit systems from interacting with one another. A transmission channel 1930 transmits a probe field represented by the coherent state |α

_(a) to the interaction region 1902. The four entangled qubit systems and the probe field input to the interaction region 1902 are represented by the input state |Ψ

_(in). The input state |ψ₄

represents four entangled two-state, matter-based qubit systems that are each represented by the qubits |ψ

₁, |ψ

₂, |ψ

₃, and |ψ

₄. The probe field separately interacts with each of the qubit systems at each of the interaction mediums 1920-1923, respectively. The probe field input to the interaction region 1902 is initially thought of as having a phase angle equal to “0” as shown in FIG. 12.

The interaction mediums 1920-1923 separately couple the probe field to each qubit system and produce an electric, or magnetic, dipole interaction between the probe field and each of the qubit systems at each of the interaction mediums. The dipole interaction at each interaction medium is characterized by an interaction Hamiltonian: Ĥ _(a,j) ^((γ))=

χ{circumflex over (σ)}_(z,j)(â ^(†) e ^(iγ) +âe ^(−iγ)) where

â^(†)e^(iγ)+âe^(−iγ) is a quadrature operator that operates on the probe field |α

_(a); and

γ is a dipole interaction parameter.

When the dipole interaction parameter γ equals “−π/2,” the quadrature operator is a Y-quadrature operator, and when the parameter γ equals “0,” the quadrature operator is an X-quadrature operator.

The dipole interaction at each of the interaction mediums 1920-1923 occurs for a period of time t and, in general, is mathematically represented by a dipole time-evolution operator: $\begin{matrix} {{\hat{U}}_{j}^{(\gamma)} = {\exp\left( {{- {\mathbb{i}}}\quad{\hat{H}}_{a,j}^{\gamma}{t/\hslash}} \right)}} \\ {= {\hat{D}\left( {{\hat{\sigma}}_{z,j}{\eta(\gamma)}} \right)}} \end{matrix}$ where

{circumflex over (D)}({circumflex over (σ)}_(z,j)η(γ))=exp({circumflex over (σ)}_(z,j)η(γ)â^(†)−{circumflex over (σ)}_(z,j)η*/(γ)â) is also called the “displacement operator;” and

η(γ)=θe^(i(γ−π/2)).

In a third embodiment of the present invention, the dipole interaction parameter γ associated with each of the interaction mediums 1920-1923 equals “−π/2.” As a result, the dipole-time evolution operator is mathematically represented by the displacement operator: ${\hat{U}}_{j}^{({\pi/2})} = {{\hat{D}\left( {{\hat{\sigma}}_{z,j}{\eta\left( {- \frac{\pi}{2}} \right)}} \right)} = {\exp\left( {\theta\quad{{\hat{\sigma}}_{z,j}\left( {\hat{a} - {\hat{a}}^{\dagger}} \right)}} \right)}}$ The qubit systems pass unchanged through each of the interaction mediums 1920-1923. The dipole interaction at the interaction mediums changes the average number of photons associated with the probe field depending on whether the qubit system transmitted to the interaction medium is in the state |0

or in the state |1

. For example, at the interaction medium 1920, a dipole interaction between a qubit system in the state |0

₁ and the probe field in the state |α

_(a) produces a displaced probe field state |α+θ

_(a), which is mathematically determined by applying the displacement operator Û₁ ^((π/2)) as follows: Û₁ ^((π/2))|0

₁|α

_(a)=|0

₁|α+θ

_(a)

The dipole interaction between the qubit system in the state |1

₁ and the probe field in the state |α

_(a) at the interaction medium 1120 produces a displaced probe field in the state |α−θ

_(a), which is also mathematically determined by: Û ₁ ^((π/2))|1

₁|α

_(a)=|1

₁|α−θ

_(a) Note that the displacement operator Û₁ ^((π/2)) does not change the state of the qubit basis states.

The interaction region 1902 receives a quantum system represented by the input state |Ψ

_(in) and outputs an entangled probe-qubit system represented by an entangled output state |Ψ

_(out). The output state |Ψ

_(out) is determined by successively applying the displacement operator associated with each interaction medium as follows. First, the dipole interaction at the interaction medium 1920 produces a first entangled probe-qubit system whose state is represented by: $\begin{matrix} {{{\hat{U}}_{1}^{({\pi/2})}\left. \Psi \right\rangle_{i\quad n}} = {\frac{1}{\sqrt{2}}\begin{pmatrix} {{\left( {{\hat{U}}_{1}^{({\pi/2})}\left. 0 \right\rangle_{1}\left. \alpha \right\rangle_{a}} \right)\left. 0 \right\rangle_{2}\left. 0 \right\rangle_{3}\left. 0 \right\rangle_{4}} +} \\ {{{\mathbb{e}}^{\mathbb{i}4\varphi}\left( {{\hat{U}}_{1}^{({\pi/2})}\left. 1 \right\rangle_{1}\left. \alpha \right\rangle_{a}} \right)}\left. 1 \right\rangle_{2}\left. 1 \right\rangle_{3}\left. 1 \right\rangle_{4}} \end{pmatrix}}} \\ {= {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha + \theta} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha - \theta} \right\rangle_{a}}} \right)}} \end{matrix}$ The first entangled probe-qubit system is input to the interaction medium 1921, which produces a second entangled probe-qubit system whose state is represented by: ${{\hat{U}}_{2}^{({\pi/2})}{\hat{U}}_{1}^{({\pi/2})}\left. \Psi \right\rangle_{i\quad n}} = {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha + {2\theta}} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha - {2\theta}} \right\rangle_{a}}} \right)}$ The second entangled probe-qubit system is input to the interaction medium 1922, which produces a third entangled probe-qubit system whose state is represented by: ${{\hat{U}}_{3}^{({\pi/2})}{\hat{U}}_{2}^{({\pi/2})}{\hat{U}}_{1}^{({\pi/2})}\left. \Psi \right\rangle_{i\quad n}} = {\frac{1}{\sqrt{2}}\begin{pmatrix} {{\left. 0000 \right\rangle\left. {\alpha + {3\theta}} \right\rangle_{a}} +} \\ {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha - {3\theta}} \right\rangle_{a}} \end{pmatrix}}$ Finally, the third entangled probe-qubit system is input to the interaction medium 1923, which produces a fourth entangled probe-qubit system whose state is represented by: $\begin{matrix} {\left. \Psi \right\rangle_{out} = {{\hat{U}}^{({\pi/2})}{\hat{U}}_{3}^{({\pi/2})}{\hat{U}}_{2}^{({\pi/2})}{\hat{U}}_{1}^{({\pi/2})}\left. \Psi \right\rangle_{i\quad n}}} \\ {= {\frac{1}{\sqrt{2}}\begin{pmatrix} {{\left. 0000 \right\rangle\left. {\alpha + {4\theta}} \right\rangle_{a}} +} \\ {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha - {4\theta}} \right\rangle_{a}} \end{pmatrix}}} \end{matrix}$

FIG. 20 shows four successive and separate phase-space diagrams 2002-2005 associated with each increase in the average number of photons acquired by the probe field at each of the interaction mediums 1920-1923, in FIG. 19, that represent an embodiment of the present invention. Each phase-space diagram in FIG. 20 includes a dashed-line disk, such as disk 2006 in diagram 2002, that represents the uncertainty disk 1208 associated with the probe field initially in the state |α

_(a), in FIG. 12. Diagrams 2002-2005 show successive increases in the average number of photons in the probe field as a result of each dipole interaction between the probe field and the qubit systems in the tensor product state 10000) at the interaction mediums 1920-1923, respectively, in FIG. 19. For example, the interaction medium 1920 receives the probe field in the state |α

_(a) and outputs the prove field in the state |α+θ

_(a) as shown in diagram 2002.

FIG. 21 shows four successive and separate phase-space diagrams 2102-2105 associated with each decrease in the average number of photons acquired by the probe field at each interaction medium, in FIG. 19, that represent embodiments of the present invention. Each phase-space diagram includes a dashed-line disk, such as disk 2106 in diagram 2102, that represents the uncertainty disk 1208 associated with the probe field, in FIG. 12. Diagrams 2102-2105 shows successive phase changes in the probe field as a result of each dipole interaction between the probe field and the qubits in the tensor product state |1111

at the interaction mediums 1920-1923, respectively, in FIG. 19. For example, the interaction medium 2020 receives the probe field in the state |α

_(a) and outputs the prove field in the state |α−θ

_(a), as shown in diagram 2021.

In FIG. 19, the interaction region 1902 outputs an entangled probe-qubit system represented by the output state: $\left. \Psi \right\rangle_{out} = {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha + {4\theta}} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha - {4\theta}} \right\rangle_{a}}} \right)}$

The Hadamard gates 1908-1911, with subsequent qubit measurements, disentangle the displaced probe fields |α+4θ

_(a) and |α−4θ

_(a) from the entangled qubit systems represented by the output state |Ψ

_(out). A Hadamard output state representing a quantum system output from the Hadamard gates 1908-1911 is mathematically determined as described above with reference to FIG. 11 and is given by: $\begin{matrix} {\left. \Psi \right\rangle_{Had} = {{\hat{H}}_{Had}^{\otimes 4}\left. \Psi \right\rangle_{out}}} \\ {= {{\left( {{\frac{1}{4}\left. \Gamma \right\rangle_{1}} + {\frac{\sqrt{6}}{4}\left. \Gamma \right\rangle_{3}} + {\frac{1}{4}\left. \Gamma \right\rangle_{5}}} \right)\left. \zeta \right\rangle_{field}^{{\pi/2} -}} + {\begin{pmatrix} {{\frac{1}{2}\left. \Gamma \right\rangle_{2}} +} \\ {\frac{1}{2}\left. \Gamma \right\rangle_{4}} \end{pmatrix}\left. \zeta \right\rangle_{field}^{{\pi/2} -}}}} \end{matrix}$ where ${\left. \zeta \right\rangle_{field}^{{\pi/2} +} = {\frac{1}{\sqrt{2}}\left( {\left. {\alpha + {4\theta}} \right\rangle + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {\alpha - {4\theta}} \right\rangle}} \right)}},{and}$ $\left. \zeta \right\rangle_{field}^{{\pi/2} -} = {\frac{1}{\sqrt{2}}{\left( {\left. {\alpha + {4\theta}} \right\rangle - {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {\alpha - {4\theta}} \right\rangle}} \right).}}$ Both of the field states |θ

_(field) ^(π/2+) and |θ

_(field) ^(π/2−) have acquired the phase shift φ and are coherent linear superposition of the probe field states output from the interaction region 1902.

The displacement operation 1904 is applied to the probe field in order to remove the average number of photons corresponding to the initial state |α

_(a) of the probe field. The displacement operation 1904 is mathematically represented by the displacement operator: {circumflex over (D)}(−α)={circumflex over (D)}*(α)=exp(−α{circumflex over (a)}+α*{circumflex over (a)}) The field states of the probe field output from the displacement operation 1904 are determined by applying the displacement operator {circumflex over (D)}(−α) to the field states and |ζ

_(field) ⁺ to give: ${{\hat{D}\left( {- \alpha} \right)}\left. \zeta \right\rangle_{field}^{{\kappa/2} +}\frac{1}{\sqrt{2}}\left( {\left. {4\quad\theta} \right\rangle + {{\mathbb{e}}^{{\mathbb{i}4}\quad\varphi}\left. {{- 4}\quad\theta}\quad \right\rangle}} \right)},{{{and}{\hat{D}\left( {- \alpha} \right)}\left. \zeta \right\rangle_{field}^{{\kappa/2} -}} = {\frac{1}{\sqrt{2}}\left( {\left. {4\quad\theta} \right\rangle - {{\mathbb{e}}^{{\mathbb{i}4}\quad\varphi}\left. {{- 4}\quad\theta} \right\rangle}} \right)}},$ respectively.

FIG. 22 is a phase-space diagram of the field states output from the displacement operation 1904, in FIG. 19, that represents an embodiment of the present invention. The displaced probe field, in diagram 2005 of FIG. 20, is translated to the probe field in the state |4θ

_(a), which is represented by directional arrow 2202 and uncertainty disk 2204. The displaced probe field, in diagram 2105 of FIG. 21, is translated to the probe field in the state |−4θ

_(a), which is represented by directional arrow 2206 and uncertainty disk 2208.

Qubit measurements performed at the qubit detectors 1912-1915 project the state of the probe field onto either the field state {circumflex over (D)}(−α)|ζ

_(field) ^(π/2+) or onto the field state {circumflex over (D)}(−α)|ζ

_(field) ^(π/2−), as described above with reference to the qubit detectors 1112-1115, in FIG. 11. In order to determine the phase shift φ acquired by the state of the probe field, a large number of quantum systems prepared in the input states |Ψ

_(in) are successively input to the implementation 1900. After each of the qubit measurements, when the probe field is in the state {circumflex over (D)}(−α)|ζ

_(field) ^(π/2+), data collected from the photon number counter 1906 are stored in a first data set, and when the probe field state is in the state {circumflex over (D)}(−α)|ζ

_(field) ^(π/2−), data collected from the photon number counter 1906 are stored in a second data set. The first data sets is used to produce a first frequency distribution associated with the field state {circumflex over (D)}(−α)|ζ

_(field) ^(π/2+), and the second data set is used to produce a second frequency distribution associated with the field state {circumflex over (D)}(−α)|ζ

_(field) ^(π/2−) as described above with reference to FIGS. 15A-15B. Note that the photon number counter 1904 absorbs the photons, and outputs a signal that can be processed by digital circuits. Photomultiplier tubes, avalanche photodiodes, and high efficiency photon detection devices, such as a visible light photon counter, can be used to determine the number of photons in the probe field.

Either the first or the second frequency distributions can be used to determine the phase shift φ. Consider, for example, using the first frequency distribution. The probability distributions associated with the probe field in the field state {circumflex over (D)}(−α)|ζ

_(field) ^(π/2+) is: $\begin{matrix} {{P^{{\kappa/2} +}(n)} = {\left\langle {n{{\hat{D}\left( {- \alpha} \right)}}\zeta} \right\rangle_{field}^{{\kappa/2} +}}^{2}} \\ {= {\frac{1}{2}{\exp\left( {{- {{4\quad\theta}}^{2}}/2} \right)}\frac{\left( {4\theta} \right)^{n}}{\sqrt{n!}}\left( {2 + {\left( {- 1} \right)^{n}{\cos\left( {4\quad\varphi} \right)}}} \right)}} \end{matrix}$ where n is the photon number. An optimal phase shift φ value can be determined by applying a least squares method, or other curve fitting method, to fit the probability distribution P^(π/2+) (n) to the data points of the first frequency distribution as a function of the phase shift φ, as described above with reference to FIG. 16. A probability distribution P^(π/2−) (n) associated with the second frequency distribution can also be used to determine the optimal phase shift φ.

In a fourth embodiment of the present invention, the parameter γ associated with each of the interaction mediums 1920-1923 equals “0.” As a result, the dipole-time evolution operator is mathematically represented by the displacement operator: Û _(j) ⁽⁰⁾⁼ {circumflex over (D)}({circumflex over (σ)}_(z,j)η(0))=exp(−θ{circumflex over (σ)}_(z,j)(â+â^(†))) The qubit systems pass unchanged through each of the interaction mediums 1920-1923. The dipole interaction at the interaction mediums changes the average number of photons associated with the probe field depending on whether the qubit system transmitted to the interaction medium is in the state |0

or in the state |1

. For example, at the interaction medium 1920, a dipole interaction between a qubit system in the state |0

₁, and the probe field in the state |α

_(a) produces a displaced probe field in the state |α+iθ)_(a), which is mathematically determined by applying the displacement operator Û₁ ⁰ as follows: Û ₁ ⁽⁰⁾|0

₁|α

_(a)=|0

₁ |α+iθ _(a) The dipole interaction between the qubit system in the state |1

₁ and the probe field in the state |α

_(a), at the interaction medium 1120 produces a displaced probe field in the state |α−iθ

_(a), which is also mathematically determined by: Û ₁ ⁽⁰⁾|1

₁|α

_(a)=|1

₁ |α−iθ _(a) Note that the displacement operator Û_(j) ⁰ does not change the state of the qubit basis states.

The interaction region 1902 receives a quantum system represented by the input state |Ψ

_(in) and outputs an entangled probe-qubit system represented by an entangled output state: $\left. \Psi \right\rangle_{out} = {\frac{1}{\sqrt{2}}\left( {{\left. 0000 \right\rangle\left. {\alpha + {{\mathbb{i}4}\quad\theta}} \right\rangle_{a}} + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 1111 \right\rangle\left. {\alpha - {{\mathbb{i}4}\quad\theta}} \right\rangle_{a}}} \right)}$

The output state |Ψ

_(out) is determined by successively applying the displacement operator associated with each interaction medium as described above.

FIG. 23 is a phase-space diagram showing the resulting displaced probe fields output from the interaction region 1902, in FIG. 19, that represent an embodiment of the present invention. Uncertainty disk 2302 and dashed-line directional arrow 2304 represent the displace probe field in the state |αa+i4θ

, and uncertainty disk 2306 and dashed-line directional arrow 2308 represent the displace probe field in the state |α−i4θ

.

In FIG. 19, the Hadamard gates 1908-1911, with subsequent qubit measurements, disentangle the displaced probe field from the entangled qubit systems represented by the output state |Ψ

_(out). A Hadamard output state is mathematically determined as described above with reference to FIG. 11 and is given by: ${{\begin{matrix} {\left. \Psi \right\rangle_{Had} = {{\hat{H}}_{Had}^{\otimes 4}\left. \Psi \right\rangle_{out}}} \\ {= {{\left( {{\frac{1}{4}\left. \Gamma \right\rangle_{1}} + {\frac{\sqrt{6}}{4}\left. \Gamma \right\rangle_{3}} + {\frac{1}{4}\left. \Gamma \right\rangle_{5}}} \right)\left. \zeta \right\rangle_{field}^{0 +}} +}} \\ {\left( {{\frac{1}{2}\left. \Gamma \right\rangle_{2}} + {\frac{1}{2}\left. \Gamma \right\rangle_{4}}} \right)\left. \zeta \right\rangle_{field}^{0 -}} \end{matrix}{where}\left. \zeta \right\rangle_{field}^{0 +}} = {\frac{1}{\sqrt{2}}\left( {\left. {\alpha + {{\mathbb{i}4}\quad\theta}} \right\rangle + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {\alpha - {{\mathbb{i}4}\quad\theta}} \right\rangle}} \right)}},{{{and}\left. \zeta \right\rangle_{field}^{0 -}} = {\frac{1}{\sqrt{2}}{\left( {\left. {\alpha + {{\mathbb{i}4}\quad\theta}} \right\rangle - {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {\alpha - {{\mathbb{i}4}\quad\theta}} \right\rangle}} \right).}}}$ The probe field, represented by the field states |ζ

_(field) ⁰⁺ and |ζ

_(field) ⁰⁻, has acquired the phase shift φ.

In FIG. 19, the displacement operation 1904 is applied to the probe field in order to remove the average number of photons associated with the initial state |α

_(a) of the probe field. The field states of the probe field output from the displacement operation 1904 are determined by applying the displacement operator {circumflex over (D)}′ (−α) to the field states |ζ

_(field) ⁰⁺ and |ζ

_(field) ⁰⁻ gives: ${{{\hat{D}\left( {- \alpha} \right)}\left. \zeta \right\rangle_{field}^{0 +}} = {\frac{1}{\sqrt{2}}\left( {\left. {\mathbb{i}4\theta} \right\rangle + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {- {\mathbb{i}4\theta}} \right\rangle}} \right)}},{and}$ ${{{\hat{D}\left( {- \alpha} \right)}\left. \zeta \right\rangle_{field}^{0 +}} = {\frac{1}{\sqrt{2}}\left( {\left. {\mathbb{i}4\theta} \right\rangle + {{\mathbb{e}}^{\mathbb{i}4\varphi}\left. {- {\mathbb{i}4\theta}} \right\rangle}} \right)}},$ respectively.

FIG. 24 is a phase-space diagram of the field states output from the displacement operation 1904, in FIG. 19, that represents an embodiment of the present invention. The displaced probe field, in FIG. 23, is translated to the probe field represented by the state |i4θ

_(a), which is represented by uncertainty disk 2402 and directional arrow 2404. The displaced probe field, in FIG. 23, is translated to the probe field represented by the state |−i4θ

_(a), which is represented by uncertainty disk 2406 and directional arrow 2408.

Qubit measurements performed at the qubit detectors 1912-1915 project the probe field output from the displacement operation 1904 onto either the field state {circumflex over (D)}(−α)|ζ

_(field) ⁰⁺ or onto the field state {circumflex over (D)}(−α)|ζ

_(field) ⁰⁻, as described above with reference to the qubit detectors 1112-1115, in FIG. 11. In order to determine the phase shift φ acquired by the probe field in the field state {circumflex over (D)}(−α)|ζ

_(field) ⁰⁺ or in the field state {circumflex over (D)}(−α)|ζ

_(field) ⁰⁻, a large number of identically prepared input states |Ψ

_(in) are successively input to the implementation 1900. After each of the qubit measurements, when the probe field is in the field state {circumflex over (D)}(−α)|θ_(field) ⁰⁺, data collected from the photon number counter 1906 is stored in a first data set, and when the probe field is in the field state {circumflex over (D)}(−α)|ζ

_(field) ⁰⁻, data collected from the photon number counter 1906 is stored in a second data set. The first data sets is used to produce a first frequency distribution associated with the field state {circumflex over (D)}(−α)|ζ^(field) ⁰⁺, and the second data set is used to produce a second frequency distribution associated with the field state {circumflex over (D)}(−α)|ζ

_(field) ⁰⁻, as described above with reference to FIGS. 15A-15B.

The first frequency distribution can be used to determine the phase shift φ by fitting a corresponding probability distribution to the first frequency distribution data set. The probability distributions associated with the field state {circumflex over (D)}(−α)|ζ

_(field) ⁰⁺ is: $\begin{matrix} {{P^{0 +}(n)} = {\left\langle {n{{\hat{D}\left( {- \alpha} \right)}}\zeta} \right\rangle_{field}^{0 +}}^{2}} \\ {= {\frac{1}{2}{\exp\left( {{- {{4\theta}}^{2}}/2} \right)}\frac{\left( {4\theta} \right)^{n}}{\sqrt{n!}}\left( {2 + {\left( {- 1} \right)^{n}{\cos\left( {4\varphi} \right)}}} \right)}} \end{matrix}$ where n is the number of photons. An optimal phase shift φ value can be determined by applying a least squares method, or other curve fitting method, to fit the probability distribution P⁰⁺ (n) to the first frequency distribution as a function of the phase shift φ as described above with reference to FIG. 16. A probability distribution P⁰⁻ (n) associated with the second frequency distribution can also be used to determine the optical phase shift φ.

In an alternate embodiment of the present invention, the phase shift φ can be extracted from a qubit system that interferes with the probe field after the Hadamard gates. FIG. 25 illustrates a fifth implementation 2500 for determining the phase shift φ acquired by the state |Ψ

₄₀₀₄ that represents an embodiment of the present invention. The implementation 2500 comprises the components described above with reference to FIG. 11 plus an interaction medium 2502, phase shift operation 2504, and qubit measurement 2506. The interaction medium 2502 is identical to the interaction mediums 1120-1123 in the interaction region 1102. At the interaction medium 2502, a qubit system in the qubit (|0

+|1

)/√{square root over (2)}, interacts with the probe field in the field states |ξ

_(field) ¹⁺ and |ξ

_(field) ¹⁻ to output an entangled probe-qubit system represented by the state: $\begin{matrix} {\left. {{\left. {\left. \Xi \right\rangle_{int} = {\frac{1}{2}\left( {{\left. 0 \right\rangle\left. {\alpha\quad{\mathbb{e}}^{\mathbb{i}5\theta}} \right\rangle_{a}} \pm {\mathbb{e}}^{\mathbb{i}4\varphi}} \right.1}} \right\rangle }\alpha\quad{\mathbb{e}}^{- {\mathbb{i}5\theta}}} \right\rangle_{a} \pm} \\ \left. \left. {{\left. {{{{{\mathbb{e}}^{\mathbb{i}4\varphi}\left. 0 \right\rangle\left. {\alpha\quad{\mathbb{e}}^{- {\mathbb{i}3\theta}}} \right\rangle_{a}} +}}1} \right\rangle }\alpha\quad{\mathbb{e}}^{\mathbb{i}3\theta}} \right\rangle_{a} \right) \end{matrix}$ The homodyne detection 1104 projects the entangled state of the entangled probe-qubit system onto a state that contains the phase shift φ. Note that the projection is an entirely random process, but the output state is determined by the homodyne detection 1104. For example, when the homodyne detection 1104 reveals the states |αe^(i5θ)

_(a) and |αe^(−i5θ)

_(a), the output state is: $\left. \left. {\left. \Xi \right\rangle_{int} = {\frac{1}{2}\left( {{{\mathbb{e}}^{{\mathbb{i}}\quad{\phi{(x)}}}\left. 0 \right\rangle} \pm {{\mathbb{e}}^{\mathbb{i}4\phi}{\mathbb{e}}^{- {{\mathbb{i}\phi}{(x)}}}}} \right.1}} \right\rangle \right),$ or when the homodyne detection 1104 reveals the states |αae^(i3θ)

_(a) and |αe^(−i3θ)

_(a), the output state is: $\left. \left. {\left. \Xi \right\rangle_{int} = {\frac{1}{2}\left( {{{\pm {\mathbb{e}}^{\mathbb{i}4\varphi}}{\mathbb{e}}^{- {{\mathbb{i}\phi}{(x)}}}\left. 0 \right\rangle} + {\mathbb{e}}^{{\mathbb{i}\phi}{(x)}}} \right.1}} \right\rangle \right)$

In either case, the phase φ(x) is identified by the homodyne detection 1104 and removed by the phase shift operation 2504. The phase shift operation 2504 is mathematically represented by the operator {circumflex over (Φ)}=exp (−iφ(x)σ_(z)/2), which produces the states: $\left. \left. {{\left. \left. {{\hat{\Phi}\quad\left. \Xi \right\rangle_{int}} = {\frac{1}{2}\left( {\left. 0 \right\rangle \pm {\mathbb{e}}^{\mathbb{i}4\varphi}} \right.1}} \right\rangle \right),{and}}{{\hat{\Phi}\quad\left. \Xi \right\rangle_{int}} = {\frac{1}{2}\left( {{{\pm {\mathbb{e}}^{\mathbb{i}4\varphi}}\left. 0 \right\rangle} +} \right.1}}} \right\rangle \right)$ The phase shift φ can then be determined via a standard single qubit interference experiment at the qubit measurement 2508.

In an alternate embodiment of the present invention, the phase shift φ can be extracted from a qubit system that interferes with the probe field after the Hadamard gates. FIG. 26 illustrates a sixth implementation 2600 for determining the phase shift φ acquired by the state |ψ₄

that represents an embodiment of the present invention. The implementation 2600 comprises the components described above with reference to FIG. 19 plus an interaction medium 2602, phase shift operation 2604, and qubit measurement 2606. The interaction medium 2602 is identical to the interaction mediums 1920-1923 in the interaction region 1902. At the interaction medium 2602, a qubit system in the qubit (|0

+|1

)/√{square root over (2)}, interacts with the probe field in the field states |ξ

_(field) ¹⁺ and |ξ

_(field) ¹⁻ to output an entangled probe-qubit system represented by the state: $\begin{matrix} {\left. \Xi \right\rangle_{int} = {\frac{1}{2}\left( {{{{\left. 0 \right\rangle\left. {\alpha + {5\theta}} \right\rangle_{a}} \pm {{\mathbb{e}}^{i\quad 4\varphi}\left. 1 \right\rangle\left. {\alpha - {5\theta}} \right\rangle_{a}}} \pm {{\mathbb{e}}^{i\quad 4\varphi}\left. 0 \right\rangle\left. {\alpha - {3\theta}} \right\rangle_{a}}} +} \right.}} \\ \left. {\left. 1 \right\rangle\left. {\alpha + {3\theta}} \right\rangle_{a}} \right) \end{matrix}$ The photon number counter 1906 projects the entangled state onto a state that contains the phase shift φ. Note that the projection is also an entirely random process, but the output state is determined by the photon number counter 1906. For example, when the photon number counter reveals an intensity of N˜25θ², the output state is: ${\left. \Xi \right\rangle_{int} = {\frac{1}{2}\left( {\left. 0 \right\rangle \pm {{{\mathbb{e}}^{i\quad 4\varphi}\left( {- 1} \right)}^{N}\left. 1 \right\rangle}} \right)}},$ or when the photon number counter 1906 reveals an intensity of N˜9θ², the output state is: $\left. \Xi \right\rangle_{int} = {\frac{1}{2}\left( {{\left( {- 1} \right)^{N}\left. 0 \right\rangle} \pm {{\mathbb{e}}^{i\quad 4\varphi}\left. 1 \right\rangle}} \right)}$ In either case, the phase factor (−1)^(N) is identified by the photon number counter 1906 and removed by the phase shift operation 2604. The phase shift operation 2604 is mathematically represented by the operator {circumflex over (Ω)}=exp(−iNπ/2)exp(iNπ{circumflex over (σ)}_(z)/2), which produces the states: ${{\hat{\Omega}\left. \Xi \right\rangle_{int}} = {\frac{1}{2}\left( {\left. 0 \right\rangle \pm {{\mathbb{e}}^{i\quad 4\varphi}\left. 1 \right\rangle}} \right)}},{and}$ ${\hat{\Omega}\left. \Xi \right\rangle_{int}} = {\frac{1}{2}\left( {{{\pm {\mathbb{e}}^{i\quad 4\varphi}}\left. 0 \right\rangle} + \left. 1 \right\rangle} \right)}$ The phase shift φ can then be determined via a standard single qubit interference experiment at the qubit measurement 2608.

Although the present invention has been described in terms of particular embodiments, it is not intended that the invention be limited to these embodiments. Modifications within the spirit of the invention will be apparent to those skilled in the art. For example, in an alternate embodiment of the present invention, rather than employing interaction regions having two or more interaction mediums as described above with reference to FIGS. 11, 17, and 19, the interaction region can comprise a single interaction medium. For example, FIG. 27 illustrates an interaction region comprising a single interaction medium that represents one of many embodiments of the present invention. In FIG. 27, the probe field is repeatedly transmitted into an interaction medium 2702 via a loop transmission channel 2704. Transmission of the probe field into the interaction medium 2702 is synchronized with the transmission of the qubit systems so that the probe interacts separately with each qubit system. In an alternate embodiment of the present invention, rather than performing the displacement operation 1904 on the probe field after the interaction region 1902, in FIG. 19, the displacement operation 1904 can be performed before the probe field is input to the interaction region 1902. In an alternate embodiment of the present invention, the displacement operation 1904 can be eliminated by using the probe field in the vacuum state |0

_(a). In an alternate embodiment of the present invention, the use of other interaction mediums characterized by different interaction Hamiltonians can be used. In an alternate embodiment of the present invention, the Hadamard gates can be replaced by x-Pauli and y-Pauli gates. In alternate embodiments of the present invention, multiple probe beams can be used to interact with subsets of qubits and then interfere at a beam-splitter prior to a measurement. In alternate embodiments of the present invention, the entangled qubit systems can be quantum dots, magnetic moments, such as superconducting flux or magnetic spin qubits, or nanoelectromechanical systems. In alternate embodiments of the present invention, each of the N matter-based qubits can be trapped within one of the N separate interaction mediums in the interaction region. The probe field |α

_(a) is transmitted to each of the N interaction mediums to separately interact with the N qubits. In an alternate embodiment of the present invention, an interaction region comprises numerous different kinds of interaction mediums. For example, an interaction region comprises a combination of interaction mediums for coupling matter-based qubits to a probe and Kerr interaction mediums for coupling photon-based qubits to the probe.

The foregoing description, for purposes of explanation, used specific nomenclature to provide a thorough understanding of the invention. However, it will be apparent to one skilled in the art that the specific details are not required in order to practice the invention. The foregoing descriptions of specific embodiments of the present invention are presented for purposes of illustration and description. They are not intended to be exhaustive of or to limit the invention to the precise forms disclosed. Obviously, many modifications and variations are possible in view of the above teachings. The embodiments are shown and described in order to best explain the principles of the invention and its practical applications, to thereby enable others skilled in the art to best utilize the invention and various embodiments with various modifications as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the following claims and their equivalents: 

1. A method for determining a phase shift acquired by an entangled N-qubit system represented by a NOON state, the method comprising: for M trials, producing an entangled probe-qubit system by coupling a probe electromagnetic field with each qubit system in the entangled N-qubit system, transferring a phase shift from the qubit systems to the probe electromagnetic field by transforming each qubit system state in the entangled probe-qubit system into a linear superposition of qubit basis states, performing a measurement on the probe electromagnetic field in order to obtain an associated measurement result, and based on the measurement result falling within a measurement-result interval of a larger continuous spectrum of measurement results, incrementing a counter associated with the measurement-result interval; producing a frequency distribution by normalizing a set of counters, each counter associated with each measurement result interval; and determining the phase shift by fitting a probability distribution associated with the probe electromagnetic field to the frequency distribution as a function of the phase shift.
 2. The method of claim 1 further comprising projecting the probe electromagnetic field onto either a first field state or a second field state by performing qubit detection measurements on each of the qubit systems.
 3. The method of claim 1 wherein the qubit systems further comprises one of: matter-based qubit systems; photon-based qubit systems; a linear superpositions of two quantum dot states; and a linear superposition of two magnetic states.
 4. The method of claim 1 wherein the coupling further comprises dipole interactions at an interaction medium represented by an interaction Hamiltonian: Ĥ=

χâ^(†)â{circumflex over (σ)}_(z) where χ is a constant that represents the coupling strength between the probe field and qubit systems; â^(†)â is a number operator that operates on the probe; and {circumflex over (σ)}_(z)=|0

0|−|1

1|.
 5. The method of claim 1 wherein the coupling further comprises a cross-Kerr interaction at an interaction medium represented by a cross-Kerr interaction Hamiltonian: Ĥ=

χâ^(†)â{circumflex over (b)}^(†){circumflex over (b)} where â^(†)â is a number operator that operates on the probe field; and {circumflex over (b)}^(†){circumflex over (b)} is a number operator that operates on a photon basis state.
 6. The method of claim 1 wherein the coupling further comprises an interaction at an interaction medium represented by an interaction Hamiltonian: {circumflex over (H)}=

χ{circumflex over (σ)}_(z)(â ^(†) e ^(iγ) +âe ^(−iγ)) where â^(†)e^(iγ)+âe^(−iγ) is a quadrature operator of the probe and is determined by the parameter γ.
 7. The method of claim 1 wherein transferring the phase shift from the qubit systems to the probe electromagnetic field further comprises transmitting each qubit system through one of the following: a Hadamard gate; an x-Pauli gate; or a y-Pauli gate.
 8. The method of claim 1 wherein performing the measurement on the probe electromagnetic field further comprises a homodyne detection measurement.
 9. The method of claim 1 wherein performing the measurement on the probe electromagnetic field further comprises a photon number counter.
 10. The method of claim 1 wherein determining the phase shift further includes employing a least squares method to fit the probability distribution to the frequency distribution as a function of the phase shift φ.
 11. A method for transferring a phase shift acquired by an entangled N-qubit system represented by a NOON state to a probe electromagnetic field, the method comprising: producing an entangled probe-qubit system by coupling a probe electromagnetic field with each qubit system in the entangled N-qubit system; and transferring a phase shift from the qubit systems to the probe electromagnetic field by transforming each qubit system state of the entangled probe-qubit system into a qubit system represented by a linear superposition of qubit basis states.
 12. The method of claim 11 further comprising projecting the probe electromagnetic field onto either a first field state or a second field state by performing qubit detection measurements on each of the qubit systems.
 13. The method of claim 11 wherein the qubit systems further comprise one of: matter-based qubit systems; photon-based qubit systems; a linear superposition of two quantum dot states; and a linear superposition of two magnetic states.
 14. The method of claim 11 wherein the coupling further comprises dipole interactions at an interaction medium represented by an interaction Hamiltonian: Ĥ=

χâ^(†)â{circumflex over (σ)}_(z) where χ is a constant that represents the coupling strength between the probe field and qubit systems; â^(†)â is a number operator that operates on the probe; and {circumflex over (σ)}_(z)=|0

0|−|1

1|.
 15. The method of claim 11 wherein the coupling further comprises a cross-Kerr interaction at an interaction medium represented by a cross-Kerr interaction Hamiltonian: {circumflex over (H)}=

χ{circumflex over (σ)}_(z)(â ^(†) e ^(iγ) +âe ^(−γ)) where â^(†)â is a number operator that operates on the probe field; and {circumflex over (b)}^(†){circumflex over (b)} is a number operator that operates on a photon basis state.
 16. The method of claim 11 wherein the coupling further comprises an interaction at an interaction medium represented by an interaction Hamiltonian: Ĥ=

χ{circumflex over (σ)}_(z)(â ^(†) e ^(iγ) +âe ^(iγ)) where â^(†)e^(iγ)+âe^(−iγ) is a quadrature operator of the probe and is determined by the parameter γ.
 17. The method of claim 11 wherein transferring the phase shift from the qubit systems to the probe electromagnetic field further comprises transmitting each qubit system through one of the following: a Hadamard gate; an x-Pauli gate; or a y-Pauli gate.
 18. A method for measuring a phase shift acquired by a probe electromagnetic field entangled with N qubit systems, the method comprising: for M trials, performing a measurement on the probe electromagnetic field in order to obtain an associated measurement result, and based on the measurement result falling within a measurement-result interval of a larger continuous spectrum of measurement results, incrementing a counter associated with each measurement-result interval; producing a frequency distribution by normalizing a set of counters, each counter associated with each measurement result interval; and determining the phase shift by fitting a probability distribution associated with the probe electromagnetic field to the frequency distribution as a function of the phase shift.
 19. The method of claim 18 wherein performing the measurement on the probe electromagnetic field further comprises a homodyne detection measurement.
 20. The method of claim 18 wherein performing the measurement on the probe electromagnetic field further comprises a photon number counter.
 21. The method of claim 18 wherein determining the phase shift further includes employing a least squares method to fit the probability distribution to the frequency distribution as a function of the phase shift φ. 